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In our previously published papers, it was proved that the chromodynamics with massive gluons 
can well be set up on the gauge-invariance principle. The quantization of the chromodynamics 
was perfectly performed in the both of Hamiltonian and Lagrangian path-integral formalisms by 
using the Lagrangian undetermined multiplier method. In this paper, It is shown that the quantum 
theory is invariant with respect to a kind of BRST-transformations. From the BRST-invariance 
of the theory , the Ward-Takahashi identities satisfied by the generating functionals of full Green 
functions, connected Green functions and proper vertex functions are successively derived. As an 
application of the above Ward-Takahashi identities, the Ward-Takahashi identities obeyed by the 
massive gluon and ghost particle propagators and various proper vertices are derived and based on 
these identities, the propagators and vertices are perfectly renormalized. Especially, as a result of 
the renormalization, the Slavnov-Taylor identity satisfied by renormalization constants is natually 
deduced. To demonstrate the renormalizability of the theory, the one-loop renormalization of the 
theory is carried out by means of the mass-dependent momentum space subtraction scheme and the 
renormalization group approach, giving an exact one-loop effective coupling constant and one-loop 
effective gluon and quark masses which show the asymptotically free behaviors as the same as those 
given in the quantum chromodynamics with massless gluons. 

RAGS: 12.38.Aw, ll.lO.Gh, 12.38.Bx, ll.lO.Jj 


I. INTRODUCTION 

According to the conventional concept of quantum chromodynamics (QCD), in order to keep the Lagrangian to 
be gauge-invariant, the gluons must be massless. On the contrary, in many previous investigations of glueballs, an 
effective gluon mass was phenomenologically introduced so as to get reasonable theoretical results [1-4]. The gluon 
mass was supposed to be generated dynamically from the interaction with the physical vacuum of the Yang-Mills 
theory [3] or through strong gluon-binding force [5]. Apparently, these arguments would not be considered to be 
stringent and logically consistent with the concept of the ordinary QCD. In our previous papers [6-8], it was argued 
that the QCD with massive gluons, as a non-Abelian massive gauge field theory in which the masses of all gauge 
fields are the same, can actually be set up on the principle of gauge-invariance without need of introducing the Higgs 
mechanism. The essential points to achieve this conclusion are as follows, (a) The gluon fields must be viewed as 
a constrained system in the whole space of vector potentials and the Lorentz condition, as a necessary constraint, 
must be introduced from the beginning and imposed on the massive Yang-Mills Lagrangian; (b) The gauge-invariance 
of a gauge field theory should be generally examined from its action other than from the Lagrangian because the 
action is of more fundamental dynamical meaning than the Lagrangian. Particularly, for a constrained system such 
as the gluon field, the gauge-invariance should be seen from its action given in the physical subspace defined by 
the Lorentz condition because the fields exist and move only in the physical subspace; (c) In the physical subspace, 
only infinitesimal gauge transformations are possibly allowed and necessary to be considered in the examination of 
whether the theory is gauge-invariant or not; This fact was clarified originally in Ref. [9]; (d) To construct a correct 
gauge field theory, the residual gauge degrees of freedom existing in the physical subspace must be eliminated by the 
constraint condition on the gauge group. This constraint condition may be determined by requiring the action to be 
gauge-invariant. Based on these points of view, it is easy to prove that the QCD with massive gluons established in 
our previous papers [6-8] is gauge-invariant. 

In Refs. [6-8], the quantization of the QCD with massive gluons (will be called massive QCD later on) was 
Lorentz-covariantly performed in the both of Hamiltonian and Lagrangian path-integral formalisms by employing 
the Lagrange undetermined multiplier method. In this paper, it will be shown that the quantum theory has an 
important property that the effective action appearing in the generating functional of Green functions is invariant 
with respect to a kind of BRST-transformations [10]. Thus, the theory is set up from beginning to end on the basis 
of gauge-invariance principle. From the BRST-symmetry of the theory, we will derive various Ward-Takahashi (W-T) 
identities [II-I7] satisfied by the generating functionals of Green functions and proper vertex functions. These W-T 
identities are of special importance in proofs of unitarity and renormalizability of the theory [18]. Furthermore, from 
the W-T identities obeyed by the generating functionals, we will derive W-T identities satisfied by the massive gluon 
propagator, the massive ghost particle propagator, the gluon three-line proper vertex, the gluon four-line proper 
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vertex and the quark-gluon proper vertex which appear in the perturbative expansion of S-matrix elements. Based 
on these W-T identities, the aforementioned propagators and vertices will be perfectly renormalized. As a result 
of the renormalizations, the Slavnov-Taylor (S-T) identity satisfied by the renormalization constants [19, 20] will be 
derived and shown to be formally the same as that given in the QCD with massless gluons (will be called massless 
QCD hereafter). This identity is much useful for practical calculations of the renormalization by the approach of 
renormalization group equation [21-23]. It should be mentioned that in the previous literature [16,24-28], the massive 
non-Abelian gauge field theory without involving Higgs bosons in it was not considered to be renormalizable and/or 
unitary. This conclusion was drawn from the theories which were not established correctly because the unphysical 
degrees of freedom involved in the theories are not eliminated at all by introducing appropriate constraint conditions. 
In our theory, the unphysical degrees of freedom appearing in the theory, i.e., the unphysical longitudinal components 
of vector potentials for the gluon fields and the residual gauge degrees of freedom existing in the subspace defined by 
Lorentz condition are respectively eliminated by the introduced Lorentz condition and the ghost equation which acts 
as the constraint condition on the gauge group. This guarantees that the massive QCD established in our previous 
papers is not only unitary, but also renormalizable. To demonstrate further the renormalizability of the theory, in 
this paper, the one loop renormalization will specifically be carried out by means of the mass-dependent momentum 
space subtraction scheme and the renormalization group equation (RGE), giving an exact one-loop effective coupling 
constant and one-loop effective quark and gluon masses without any ambiguity. 

The arrangement of this paper is as follows. In section 2, we will derive the BRST-transformations under which the 
effective action of the massive QCD is invariant. In section 3, we will derive the W-T identities satisfied by various 
generating functionals. In section 4, to illustrate applications of the above W-T identity, the W-T identities obeyed by 
the massive gluon propagator and ghost propagator will be derived and the renormalization of these propagators will 
be discussed. In sections 5-7, the W-T identities obeyed by the gluon three-line vertex, the gluon four-line vertex and 
the quark-gluon vertex will be derived and the renormalizations of the vertices will be discussed, respectively. Section 
8 serves to derive the one-loop effective coupling constant and the effective gluon mass. Section 9 is used to derive 
the one-loop effective quark masses. In the last section, some conclusions and discussions are made. In Appendix, the 
W-T identities will be given by an alternative derivation. 


II. BRST- TRANSFORMATION 

In Ref. [7], the QCD with massive gluons is set up by starting from the Lagrangian 

L = ^{* 7 ^( 9 ^ - ( 2 . 1 ) 

where denotes the quark field function , 'tp{x) is its Dirac-conjugate, T® = A“/2 are the color matrices (the 
generators of gauge group SU (3)), m is the quark mass and M is the gluon mass. The above Lagrangian is constrained 
by the Lorentz condition 


d^Al = 0 . ( 2 . 2 ) 

Under this condition, as was proved in Refs. [6-8], the action given by the Lagrangian in Eq. (2.1) is invariant with 
respect to the following gauge transformations: 


5Al = ^DfC\ 

d-ipix) = ig^T°-C°‘(x)4’ix), (2.3) 

6'ip{x) = ig^'ip{x)T^C°'{x) 

where 

Df = 5^'^d^-gr^^Al (2.4) 

is the covariant derivative. In the above, we have set the parametric functions of the gauge group 9°'{x) = ^C“(x) in 
which ^ is an infinitesimal Grassmann number and C^^x) are the ghost field functions. 

The quantization of the massive QGD was carried out by different approaches in Ref. [6]. A simpler quantization 
is performed in the Lagrangian path-integral formalism by means of the Lagrange undetermined multiplier method 
which was shown to be equivalent to the Faddeev-Popov approach of quantization [9]. For this quantization, it is 
convenient to generalize the QCD Lagrangian and the Lorentz condition to the following forms: 
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(2.5) 


Cx = - m}4, - 

+iM 2 y 4 “^^“ - ia(A “)2 


and 


(9^.4“ + aA“ = 0 (2.6) 

where A“(x) are the extra functions which will be identified with the Lagrange multipliers and a is an arbitrary 
constant playing the role of gauge parameter. According to the general procedure for constrained systems, the 
constraint in Eq. (2.6) may be incorporated into the Lagrangian in Eq. (2.5) by the Lagrange multiplier method, 
giving a generalized Lagrangian such that 

Lx = - igT<^Al) - m}^ - 

+A“9^A“ + ia(A“)2. 

This Lagrangian is obviously not gauge-invariant. However, for building up a correct gauge field theory, it is necessary 
to require the dynamics of the system, i.e. the action given by the Lagrangian (2.7) to be invariant under the gauge 
transformations denoted in Eq. (2.3). By this requirement, noticing the identity _ q applying the 

constraint condition in Eq.(2.6), we find 

5Sx = -- f d^xaM“(x)5^(I?f (a;)C'^(x)) = 0 (2.8) 

a J ^ 


where 


I?“'(x)=5“*'^a^ + <(x) (2.9) 

*—‘ai 

in which cr^ = aM^, is the D’Alembertian operator and D^{x) was defined in Eq. (2.4). From Eq. (2.6) we see 
A9'^A“ = —A“ 0. Therefore, to ensure the action to be gauge-invariant, the following constraint condition on the 
gauge group is necessary to be required 


d>i,{Vf{x)C\x)) = 0 (2.10) 

which usually is called ghost equation. When this constraint condition is incorporated into the Lagrangian in Eq. 
(2.7) by the Lagrange multiplier method, we obtain a more generalized Lagrangian as follows 

Lx = - igT-AD - 

-hA“a^A“ -h ia(A“)2 -y ^ ’ 


where C'“(x), acting as Lagrange undetermined multipliers, are the new scalar variables conjugate to the ghost 
variables C°‘{x). 

As we learn from the Lagrange undetermined multiplier method, the dynamical and constrained variables as well 
as the Lagrange multipliers in the Lagrangian (2.11) can all be treated as free ones, varying arbitrarily. Therefore, 
we are allowed to use this kind of Lagrangian to construct the generating functional of Green functions 


77,77,^“,^“] = i/L>(A“,V>^ 7 />,C'“,C'“, A“)ea:p{i/dS[/:A(a;) 
+^^^{x)Al^{x)_+ V'?? -F ryiA + C{x)C^{x) -h C'“(x)e(a;)]} 

= ^ d'^x[J°'^^{x)Al{x) 

+ipri + ri'tl; + ^“(x)C'“(x) + (^“(a;)^®)^)]} 


where the last equality is obtained by carrying out the integral over A“(a;), H(A“,---,A“) denotes the functional 

integration measure, and are the external sources coupled to the gluon, quark and ghost fields, N is the 

normalization constant and 


S = J - igT’^A^) - m}V; - 

-f iM2A“^A“ - ^(9^A“)2 - 

is the effective action given in arbitrary gauges. 


(2.13) 


3 



Similar to the massless QCD, for the massive QCD, there are a set of BRST-transformations including the in¬ 
finitesimal gauge transformations shown in Eq. (2.3) and the transformations for the ghost fields under which the 
effective action is invariant. The transformations for the ghost fields may be found from the stationary condition of 
the effective action under the BRST-transformations. By applying the transformations in Eq. (2.3) to the action in 
Eq. (2.13), one can derive 

SS = J = 0 

This expression suggests that if we set 

5 ( 7 “ = 

a 

and 

d^^6{VfC^) = 0. (2.16) 

The action will be invariant. Eq. (2.15) gives the transformation law of the ghost field variable C°‘{x) which is the 
same as the one in the massless gauge field theory. From Eq. (2.16), we may derive a transformation law of the ghost 
field variables C°‘{x). Noticing the relation in Eq. (2.9), we can write 

5{Vf{x)C\x)) = ^dldC^^ix) + 5{Df{x)C\x)). (2.17) 

In the massless QCD, it has been proved that [13-17] 

6{Df{x)C\x)) = Df{x)[5C\x) + \gf^<^C^{x)C\x)]. (2.18) 

With this result, Eq. (2.17) can be written as 

6{v;\x)C\x)) = Vf{x)5C\x) - Df{x)5C^,{x) 

where 

5C^{x) = -^/“*'“C'^(x)C'“(a;). 

On substituting Eq. (2.19) into Eq. (2.16), we have 

M^\x)6C\x) = Mg\x)5Cl^{x) 

where we have defined 

M^\x) = d^v;\x) = + a^) - (x)a^( (2.22) 

and 

Mg\x) = d^Df{x) = M^\x) - (2.23) 

It is noted that the operator in Eq. (2.22 ) is just the operator appearing in Eq. (2.10). Corresponding to Eq.(2.10), 
we may write an equation satisfied by the Green function A“*'(a; — y) 

M‘^^{x)A^\x -y) = 5“'’5^(x - y). (2.24) 

The function A°‘^{x — y) is nothing but the exact propagator of the ghost field which is the inverse of the operator 
M°'^{x). In the light of Eq. (2.24) and noticing Eq. (2.23), we may solve out the 5C°'{x) from Eq. (2.21) 

5C“(x) = (M-iMo5Co)“(x) = {M-i(M - u2)5Co}“(x) 

= 6Cg{x) - J d‘^yA‘^'^{x-y)6C^{y). 

This just is the transformation law for the ghost field variables C“(a:). When the gluon mass M tends to zero, 
Eq. (2.25) immediately goes over to the corresponding transformation given in the massless gauge field theory. It 
is interesting that in the Landau gauge (a = 0), due to ct = 0, the above transformation also reduces to the form 
as given in the massless theory. This result is natural since in the Landau gauge, the gluon field mass term in the 
action is gauge-invariant. However, in general gauges, the mass term is no longer gauge-invariant. In this case, to 
maintain the action to be gauge-invariant, it is necessary to give the ghost field a mass a so as to counteract the 
gauge-non-invariance of the gluon field mass term. As a result, in the transformation given in Eq. (2.25) appears a 
term proportional to . 


(2.19) 

( 2 . 20 ) 
( 2 . 21 ) 


(2.14) 

(2.15) 
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III. WARD-TAKAHASHI IDENTITIES 


This section is devoted to deriving the W-T identities for massive QCD on the basis of the BRST-symmetry of the 
theory. Since the derivations are much similar to those for the QCD with massless gluons, we only need here to give a 
brief description of the derivations. When we make the BRST-transformations shown in Eqs. (2.3), (2.15) and (2.25) 
to the generating functional in Eq. (2.12) and consider the invariance of the generating functional, the action and the 
integration measure under the transformations (the invariance of the integration measure is easy to check), we obtain 
an identity such that [11-17] 


i /I?(A“,C'“,C“,V',V’) I + dC‘'{x)^^{x) + C{x)dC^{x) 

+fi{x)S'>p{x) + (5'i/'(a;)77(x)}e*'®+-®‘®^ 

= 0 


(3.1) 


where EST is an abbreviation of the external source terms appearing in Eq. (2.12). The Grassmann number ^ 
contained in the BRST-transformations in Eq. (3.1) may be eliminated by performing a partial differentiation of Eq. 
(3.1) with respect to As a result, we get a W-T identity as follows 


i / I?(A“, ^) / d^x{J‘^^^ix)AA-ix) + AC“(a;)r (cc) - ?“(a;)AC“(a:) 

—r]{x)A'tj;{x) + A'tp{x)ri{x)}e^^^^^'^ 

= 0 


(3.2) 


where 


AAlix) = Df{x)C\x), 

AC'“(x) = ^^^^Al{x), 

AC^(x) = f - y) - a^A‘‘‘’(x - y)]AC'^(y), 

AC^(y) =-^gf^<^C^(y)C%), 

A^l){x) = igT°‘C°'{x)'ip{x), 

Ai/j(x) = ig'ip{x)T°-C°‘\x). 

These functions dehned above are hnite. Each of them differs from the corresponding BRST-transformation written 
in Eqs. (2.3), (2.15) and (2.25) by an infinitesimal Grassmann parameter 

In order to represent the composite held functions AA“, AC“, A-ip and Aij) in Eq. (3.2) in terms of differentials of 
the functional Z with respect to external sources, we may, as usual, construct a generalized generating functional by 
introducing new external sources (called BRST-sources later on) into the generating functional written in Eq. (2.12), 
as shown in the following [13-16] 


= A J X)[A“, (7“, (7“, V’, ^l}]exp{iS + i f d‘‘x[u“''AA“ -h v^AC"^ (3.4) 

-t A^C + CAV' -k -h + V^ + -ipv]} 

where u“, c and c are the sources which belong to the corresponding functions AA“ AC“, A'k and Atp, 
respectively. Obviously, the and AA“ are anticommuting quantities, while, the C, AG“, Aijj and A'ljj are 

commuting ones. We may start from the above generating functional to re-derive the W-T identity. In order that the 
identity thus derived is identical to that as given in Eq. (3.2), it is necessary to require the BRST-source terms UiA^i, 
where Ui = u°‘^, oi C, and A$i = AA“, AC°‘, Ad^ or A^k to be invariant under the BRST-transformations. 

How to ensure the BRST-invariance of the source terms? For illustration, let us introduce the source terms in such a 
fashion 

f -|- v°-SC°' + (6tjj + 6tpQ 

= f d^x[u‘^/^AA^ + v“AG“ -k CAV' -k A^C] 

where 

= = k = ?e, k = -kT 

These external sources are defined by including the Grassmann number ^ and hence products of them with ^ vanish. 
This suggests that we may generally define the sources by the following condition 


(3.5) 

(3.6) 
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Uli = 0. 


(3.7) 


Considering that under the BRST-transformation, the variation of the composite field functions given in the general 
gauges can be represented in the form where d’i are functions without including the parameter clearly, 

the definition in Eq. (3.7) for the sources would guarantee the BRST- invariance of the BRST-source terms. When 
the BRST-transformations in Eqs. (2.3), (2.15) and (2.25) are made to the generating functional in Eq. (3.4), due 
to the definition in Eq. (3.7) for the sources, we have Ui5/S.^i = 0 which means that the BRST-source terms give a 
vanishing contribution to the identity in Eq. (3.1). Therefore, we still obtain the identity as shown in Eq. (3.1) except 
that the external source terms is now extended to include the BRST-external source terms. This fact indicates that 
we may directly insert the BRST-source terms into the exponent in Eq. (3.1) without changing the identity itself. 
When performing a partial differentiation of the identity with respect to we obtain a W-T identity which is the 
same as written in Eq. (3.2) except that the BRST-source terms are now included in the identity. Therefore, Eq. 
(3.2) may be expressed as 


+^(^)3cf7) 




(x) 

= 0 . 


- V{x)- 


Sv‘^{x) SC{x) 

' ■■■,(] 


(3.8) 


This is the W-T identity satisfied by the generating functional of full Green functions. 

On substituting in Eq. (3.8) the relation [13-16] 

Z = (3.9) 

where W denotes the generating functional of connected Green functions, one may obtain a W-T identity expressed 
by the functional W 


= 0 . 


(3.10) 


Erom this identity, one may get another W-T identity satisfied by the generating functional T of proper (one-particle- 
irreducible) vertex functions. The functional T is usually defined by the following Legendre transformation [13-16] 


^[7l“^ K, U] = w[j-, f, e, fj, m C, C] 

- / d‘‘a;[ -h + + ^v] 

where and ij} are the field variables defined by the following functional derivatives 

Aa{ry\ _ SW 

i’i^) = -4m. = 4m- 

Erom Eq.(3.11), it is not difficult to get the inverse transformations [13-16] 


J“^(x) = — ST ca(„\ _ 3r 


= -4r- 


■ SA^ix) ’ ^ — SC-^ix 

fj{x) = 


It is obvious that 


5i/>( 


6w _ 6r SW _ sr sw _ jt sw _ jt 

<5^“ ’ Sv°' (5v“’ i5C ’ 


Employing Eqs. (3.13) and (3.14), the W-T identity in Eq. (3.10) will be written as [13-16] 

f AAryS _§T_ST_ I <5r 3r I 3r 3r 

1 1 5A^{x) 5u^f^{x) 5C^{x) 6v^{x) S'tp{x) S(^{x) 

^ S^Pix) 5C(x) ^ sc"{x)^ 


= 0 . 


This is the W-T identity satisfied by the generating functional of proper vertex functions. 


(3.11) 


(3.12) 


(3.13) 


(3.14) 


(3.15) 
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The above identity may be represented in another form with the aid of the so-called ghost equation of motion. The 
ghost equation may easily be derived by hrstly making the translation transformation: (7“ ^ -I- A“ in Eq. (2.12) 

where A“ is an arbitrary Grassmann variable, then differentiating Eq. (2.12) with respect to the A“ and hnally setting 
A^ = 0. The result is [13-16] 

1 J = 0. (3.16) 


When we use the generating functional defined in Eq. (3.4) and notice the relation in Eq. (2.9), the above equation 
may be represented as [13-16] 




5u°-^^{x) 


— ICF 


5 S. {x) 




(3.17) 


On substituting the relation in Eq. (3.9) into the above equation, we may write a ghost equation satisHed by the 
functional W such that 


r(x) + a(( 


5W 


6u°->^{x) 


, 6W 
Wix) 


= 0 . 


(3.18) 


From this equation, the ghost equation obeyed by the functional T is easy to be derived by virtue of Eqs. (3.12) - 
(3.14) [13-16] 


(5C“(x) (5 m“^(x) 


^^^“(x) = 0. 


(3.19) 


Upon applying the above equation to the last term in Eq. (3.15). the identity in Eq. (3.15) will be rewritten as 


J ' 


l 5A<‘, 


V ot ^ J 

= 0 . 


(3.20) 


Now, let us define a new functional E in such a manner 

r = r + ^ I d^xid^A^f. 

From this definition, it follows that 


(5r 

5A- 


r-p -I 

+ -d^^d''Ai. 




When inserting Eq. (3.21) into Eq. (3.20) and considering the relation in Eq. (3.22), we arrive at 


(3.21) 


(3.22) 


J (5A“ 


,5f 5f 
(5(7“ (5w“ 


dip S( 


(5i/j S(p 


+ M'^d’^AlC^} = 0. 


The ghost equation represented through the functional f is of the same form as Eq. (3.19) 


6C°‘{x) ^ (5m“^(x) 


(t2(7“(x) = 0. 


In the Landau gauge, since ct = 0 and = 0, Eqs. (3.23) and (3.24) respectively reduce to [13-16] 


(3.23) 


(3.24) 


and 


[ rA liLJL 

J 5A<p 


(5f (5f 
(5(7“ (5x“ 


Sip (5C Sip S( 


(5(7“ Su'^f^ 


= 0 . 


(3.25) 


(3.26) 


These equations formally are the same as those for the massless QCD. 

From the W-T identities formulated above, we may derive various W-T identities obeyed by Green functions and 
vertices, as will be illustrated in the next sections. 
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IV. GLUON AND GHOST PARTICLE PROPAGATORS 


In this section, we plan to derive the W-T identities satisfied by the massive gluon and ghost particle propagators 
by starting from the W-T identity represented in Eq. (3.8) and the ghost equation shown in Eq. (3.17) and then 
discuss their renormalization. Let us perform differentiations of the identities in Eqs. (3.8) and (3.17) with respect to 
the external sources ^“(x) and ^^{y) respectively and then set all the sources except for the source J“(x) to be zero. 
In this way, we obtain the following identities 



(4.1) 

and 


S^Z[U.u] 1 , 2 1 

+S‘^^S^{x - y)Z[J] =0. 

(4.2) 

Furthermore, on differentiating Eq. (4.1) with respect to J^{y) and then letting the source 
identity which is, in operator representation, of the form [13-16] 

J vanish, we may get an 

-dli < 0+lr[i“(x)i^(y)]l0- >=< 0+\T*[C-ix)Dt‘^iy)C‘^iy)]\0- > 
a ^ 

(4.3) 

where A“(x), C°‘{x) and C'“(x) stand for the gluon field and ghost field operators and T* symbolizes the covariant 
time-ordering product. When the source J is set to vanish, Eq. (4.2) gives such an equation [13-16] 

id!^ < Q+\T*{6-{x)bl^{y)C<^{y)}\Q- > 

+ia‘^ < 0+lr[C'“(x)C''’(y)]l0- >= 5“'’(5‘‘(x - y). 

(4.4) 

Upon inserting Eq. (4.3) into Eq. (4.4), we have 


dxdyD'^’Kx -y)- (a(T^A“'’(x - y) = -aS^^6^{x - y) 

(4.5) 

where 


iDll{x-y) =< 0+lT{i“(x)i^(y)}l0- > 

(4.6) 

which is the familiar full gluon propagator and 


iA“''(x - y) =< 0+lT{C'“(x)(7''(j/)}l0- > 

(4.7) 


which is the full ghost particle propagator. Eq. (4.5) just is the W-T identity respected by the gluon propagator 
which establishes a relation between the longitudinal part of gluon propagator and the ghost particle propagator. 
Particularly, in the Landau gauge (a = 0), as we see, Eq. (4.5) reduces to the form which exhibits the transversity of 
the gluon propagator. By the Fourier transformation, Eq. (4.5) will be converted to the form given in the momentum 
space as follows 


k>^k''Df^{k) - aa^A‘^\k) = -a5^\ (4.8) 

The ghost particle propagator may be determined by the ghost equation shown in Eq. (4.4). However, we would 
rather here to derive its expression from the Dyson-Schwinger equation [29] satisfied by the propagator which may be 
established by the perturbation method. 

A^\k) = Al\k) + Ag“'(fc)H“''''(fc)A'''^(fc) (4.9) 


where 




iAg''(A:) = i(5“'’Ao(fc) 


k'^ — a'^ + is 


(4.10) 



is the free ghost particle propagator which can be derived from the generating functional in Eq. (2.12) by a perturbative 
calculation and —denotes the proper self-energy operator of ghost particle. From Eq. (4.9), it 
is easy to solve that 




__ 

k‘^[l + Cl{k'^)] — + ie 


where the self-energy has properly been expressed as 

n{k) = k^Cl(k^). 


(4.11) 


(4.12) 


Similarly, we may write a Dyson-Schwinger equation for the gluon propagator by the perturbation procedure [29] 

D^u{k) = DIM + (4.13) 

where the color indices are suppressed for simplicity and 




■sabt dpi^ - kM/k'^ , akM/k'^ i 
k"^ — + ie — cr^ -I- ie 


(4.14) 


is the free gluon propagator which can easily be derived from the perturbative expansion of the generating functional 
in Eq. (2.12) and —in“^(fc) = —i(5“^n^^(fc) stands for the gluon proper self-energy operator. Let us decompose the 
propagator and the self-energy operator into a transverse part and a longitudinal part: 


L»A‘^(fc) = L»^^(fc) -H D'^‘'{k)M‘'ik) = Mik) + Mik) (4.15) 

where 

M{k) = ■p^''(A:)nT(fc2), n^^(fc) = ■p^^(fc)nL(fc2) 


here7^y^(fc) = aadV'^{k) = are the transverse and longitudinal projectors respectively. Considering 

these decompositions and the orthogonality between the transverse and longitudinal parts, Eq. (4.13) will be split 
into two equations 


DTp,v{k) — D^^^{k) + D^^^{k)Il^^{k)DTpi^ik) 

and 

DlM) = DlM) + 

Solving the equations (4.17) and (4.18), one can get 

jjjab — „.yab r 9^^’' ~ kpku/k"^ _^_ akpku/k'^ _, 

’ ^k‘^ + nM)-M-^ + ie k’^ + anM)-a'^ + iek 


With setting 

nM) = fc2ni(fc2) + 

which follows from the Lorentz-covariance of the operator IiT{k^) and 

anM) = k'^i\M), 


(4.17) 

(4.18) 

(4.19) 

(4.20) 

(4.21) 


Eq. (4.19) will be written as 


iDfM = 


gpy - kM/k'^ 


+ 


akpki^/k'^ 


k^[l + ni(A:2)] - M2[l - n 2 (fc 2 )] + ie k^l + UM)] -a^+ie 


(4.22) 


We would like to note that the expressions given in Eqs. (4.12), (4.20) and (4.21) can be verihed by practical 
calculations and are important for renormalizations of the propagators and the gluon mass. 

Substitution of Eqs. (4.11) and (4.22) into Eq. (4.8) yields 
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(4.23) 


nL(fc2) 


fc2[l + fi(fc2)]’ 


From this relation, we see, either in the Landau gauge or in the zero-mass limit, the ni,(fc^) vanishes. 

Now let us discuss the renormalization. The function in Eq. (4.11) and the functions ni(fc2)^ n 2 (fc^) and 

IlL(fc^) in Eq. (4.22) are generally divergent in higher order perturbative calculations. According to the conventional 
procedure of renormalization, the divergences included in the functions Cl{k'^), ni(fc2)^ n 2 (fc^) and liLik^) may be 
subtracted at a renormalization point, say, k^ = y?. Thus, we can write [13-17] 


^(fc2) = i7(/x2) + ni(fc2) = Hi (^ 2 ) ^ n^(fc 2 ), 

n2(fc2) = n2(Ai2) + n^(fc2), n^ik^) = fii(/x2) + n-(fc^) 


where Hi(^2)^ n 2 (^^), nL(/x^) and Vt‘^{kF‘), IVi{k’^), li'^ik'^), nf,(fc^) are respectively the divergent parts and the 

finite parts of the functions n(fc^), ni(fc2)^ n 2 (fc^) and IlL(fc^). The divergent parts can be absorbed in the following 
renormalization constants defined by [13-17] 


4-1 = 1 + ^ 3-1 = 1 + = l + UL{^i ^), 

z-^ = = V[i-ni(/x2)][i-n2(/x2)] 


where Z 3 and Z 3 are the renormalization constants of gluon and ghost particle propagators respectively, Z^ is the 
additional renormalization constant of the longitudinal part of gluon propagator and Zm is the renormalization 
constant of gluon mass. With the above definitions of the renormalization constants, on inserting Eq. (4.24) into Eqs. 
(4.11) and (4.22) , the ghost particle propagator and gluon propagator can be renormalized, respectively, in such a 
manner 


and 


where 


and 


iA^^’ik) = ZsiA'^^ik) 


zD-Jlik) = Z^iD^-l{k) 


iAf{k) = 




A:2[l + r2fl(fc2)] _ (j2^ _l_ 


^DZAk) = 


gfiv - kf_tk„/k'^ 


Z'^aRk^k^/k"^ 


k‘^-Ml + Ulik^) + is fc2[l + n^(A:2)] -aji + is 


} 


(4.26) 


(4.27) 


(4.28) 


(4.29) 


are the renormalized propagators in which Mr, ur and ur are the renormalized masses, aR is the renormalized 
gauge parameter, ^lR{k‘^),\\^{k‘^) and n)^(fc^) denote the finite corrections coming from the loop diagrams. They are 
defined as 


^R{k^) 


Mr = Zf^M, aR = Z 3 ur = <jr = y Z^a, 

z^iCi^P), nliP) = z^ienne) + M^n-^iP)], n^(P) = z'fii(P). 


(4.30) 


The finite corrections above are zero at the renormalization point /i. As we see from Eq. (4.29), the longitudinal part 
of the gluon propagator, except for in the Landau gauge, needs to be renormalized and has an extra renormalization 
constant Zy This fact coincides with the general property of the massive vector boson propagator (see Ref. (16), 
Chap.V). From Eqs. (4.23)-(4.25) , it is easy to find that the longitudinal part in Eq. ( 4.22) can be renormalized as 


where 


a 

fc2[l -I- flR^k"^)] — + ie 


ZsaRll + Ofl(A:^)]A/{(fc^) 


AR{k^) 


1 

fc2[l -I- fl^(fc2)] — aj^ + ie 


(4.31) 


(4.32) 
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which appears in Eq. (4.28) and the renormalization constant Zg can be expressed as 


If choosing = a r, we have 


^3 


[1 + 


4 ( 1 -^ 3 ) 

Z 3 


-1 


^3 


= ^ 3 . 


(4.33) 


(4.34) 


V. GLUON THREE-LINE VERTEX 

The aim of this section is to derive the W-T identity satisfied by the gluon three-line proper vertex and discuss its 
renormalization. For this purpose, we first derive a W-T identity satisfied by the gluon three-point Green function. 
Let us begin with the derivation from the W-T identity in Eq. (4.1) and the ghost equation in Eq. (4.2). By taking 
successive differentiations of Eq. (4.1) with respect to the sources J^iy) and T^(z) and then setting the sources to 
vanish, one may obtain the W-T identity obeyed by the gluon three-point Green function which is written in the 
operator form as follows 


-d^,Gf:^{x,y,z) = < 0+\T*[C-{x)bl\y)C'^{y)Al{z)]\0- > 

+ < Q+\T*[CAx)Al{y)bf{z)C’^{z)]\0- > (5.1) 

where 

Gt:Ax,y,z) =< Q+\nAl{x)Al{y)Al{z)]\Q- > (5.2) 

is the three-point Green function mentioned above. The identity in Eq. (5.1) will be simplihed by a ghost equation 
which may be derived by differentiating Eq. (4.2) with respect to the source TJ(^) 

< 0+\T*{DlAx)C<^{x)CAy)Al{z)}\0- > 

< o+\T[C‘^{x)dAy)A‘i{z)]\o-] >= 0. 


Taking derivatives of Eq. (5.1) with respect to y and z and employing Eq. (5.3), we get 

d^dp^Gli^ix, y, z) = aa^{dlG^<^l{z, x, y) + d^G^^Hy, x, z)} 


where 


(5.4) 


G“'’;(x,y,z) =< 0+|T{G“(x)(7'>(2/)i^(z)}|0- > . (5.5) 

In the Landau gauge or the zero-mass limit (ct = 0), Eq. (5.4) reduces to 

aoa,"a,^G“tA(uy,^) = o (5.6) 

which shows the transversity of the Green function. From Eq. (5.4), we may derive a W-T identity for the gluon 
three-line vertex. For this purpose, it is necessary to use the following one-particle-irreducible decompositions of the 
Green functions which can easily be obtained by the well-known procedure [13-16] 


Gf^A^, y,z) = f Ax'Ay'Az'zD-p(x - x') 
y.iD'bUy - y'ADliiz - (x', y', z') 


and 


G‘^^l{x,y,z) = JAx'Ay'Az'iA^^’ix - 
xiA^'Ay’ -y)iDi%{z' 


x')r“'^'‘='’'^'(x',y',z') 


(5.8) 


where iD““/(x — x') and (x — x') are respectively the gluon and the ghost particle propagators discussed in the 

preceding section, V^^(x,y^z) and r“^^(x,y,z) are the three-line gluon proper vertex and the three-line ghost-gluon 
proper vertex respectively. They are defined as [13-16] 
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(5.9) 


and 




S^T 


SAl^(x)dA'i,{y)5Al{z) 


J=0 


^abc 


ix,y,z) = 




i6C°‘{x)5C^{y)5A^^{z) 


|j=o 


(5.10) 


where J stands for all the external sources. Substituting Eqs. (5.7) and (5.8) into Eq. (5.4) and transforming Eq. 
(5.4) into the momentum space, one can derive an identity which establishes the relation between the longitudinal 
part of three-line gluon vertex and the three-line ghost-gluon vertex as follows 


= -^x(p")[x(fc")g"A“^(fc,P,g) 

+xiq'^)k^A^°'l{q,p, k)] 


(5.11) 


where we have defined 


k) = (27r)454(p + q + k)Af^\{p, q, k), 
T°-\{p, q, k) = {211^5'^{p -b <7 -b fc)A“*'^(p, < 7 , k) 


and 

x(p^) = {A^[l + Al{p^)] - -b ie}{fc^[l -b fl(p^)] - -b ie}~^ 
= [l + f2(fc2)]-i 


(5.12) 


(5.13) 


here liz,(p^) and fi(p^) are the self-energies defined in Eqs. (4.12) and (4.21). The second equality is obtained by 
inserting the relation in Eq. (4.23) into the first equality. 

Obviously, in the Landau gauge, Eq. (5.11) reduces to 

p>^q''k>'Af,\{p,q,k) = Q (5.14) 

which implies that the vertex is transverse in this case. In the lowest order approximation, owing to 


X(p") = 1 (5.15) 

and 

A(0)“'’;(p,g,fc)=5r'’% (5.16) 

where are the structure constants of the gauge group, the right hand side (RHS) of Eq. (5.11) vanishes, therefore, 
we have 

pt^q-k^A^°^ll\ip,q,k) = 0. (5.17) 

This result is consistent with that for the bare three-line gluon vertex given by the Feynman rule. 

Now, let us discuss renormalization of the three-line gluon vertex. From the renormalization of the gluon and ghost 
particle propagators described in Eqs. (4.26) and (4.27) and the definitions of the propagators written in Eqs. (4.6) 
and (4.7), one can see 


A^^ix) = V^A“/(x), 

C^ix) = V^C'^(x), C^ix) = V^Q(x) ^ ^ 

(hereafter the subscript R marks renormalized quantities). According to above relations and the definitions given in 
Eqs. (5.9), (5.10) and (5.12), we find 


K'IMp, fo k) = q, k), 

A-\ip, q, k) = q, k). 

Applying these relations, the renormalized version of the identity written in Eq. (5.11) will be 
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(5.20) 


p^q-'k^AR;tl{p,q,k) = -^XRiP^){XRik^)q’^A^<^-^{k,p,q) 
where an and an were defined in Eq. (4.30) and 


XRik^) 


1 

l + nnik^)] 


(5.21) 


is the renormalized expression of the function x{k‘^)- Iii f^e above, we have considered 


X(fc2) = Z^XRik^) (5.22) 

which follows from f)(fc^) = 0(/x^)+ f2'^(A:^), l + and nRlk"^) = Zz^^^k"^) defined in the preceding section. 

At the renormalization point chosen to be = k^ = p?, we see, Xr{p^) = 1- this case, the renormalized 

ghost-gluon vertex takes the form of the bare vertex so that the RHS of Eq. (5.20) vanishes, therefore, we have 

p^g'^fc^A^“()')^(p,(7,A:)|p2^,2^fc2^^2 =0. (5.23) 

Ordinarily, one is interested in discussing the renormalization of such three-line vertices that they are defined from 
the vertices defined in Eqs. (5.9) and (5.10) by extracting a coupling constant g. These vertices are denoted by 
A“()‘y(p, q, k) and q, k). Commonly, they are renormalized in such a fashion [13-17] 


AfA{p, q, k) = Z^ k), 

A-b^^{p,q,k) = Z^^A]^\{p,q,k) 


where Zi and Zi are referred to as the renormalization constants for the gluon three-line vertex and the ghost-gluon 
vertex, respectively. It is clear that the W-T identity shown in Eq. (5.11) also holds for the vertices A“^'y(p, g, fc) 
and A“^^(p, q, k). So, when the vertices A“^'y (p, q, k) and A“^^(p, q, k) in Eqs. (5.11) are replaced by A“^'y(p, q, k) and 

A°‘\{p,q,k) respectively and then Eq. (5.24) is inserted to such an identity, we obtain a renormalized version of the 
identity as follows 


pf^q''k^Aj^‘^ll{p,q,k) = Jxfl(p^)[XK(fc^) (5 25 ) 

x<7‘^A«“^(fc,p,g) +Xfl(9")fc^A^^=(g,p,fc)]. 

When multiplying the both sides of Eq. (5.25) with a renormalized coupling constant gn and absorbing it in the 
vertices, noticing 


AR;i^x {p, q, k) = gnAt^vx (P^ k ), 

At^'xip^ q^ k) = gRA-t^'xip^ q^ k), 

we have 

p>^q’'k^Aji^l^^{p,q,k) = -^^^XR{p'^)[XR{k^) )5 27 ) 

xq''Aj,--^{k,p,q) + XR{qnk^A^»-^{q,p,k)]. 

In comparison of Eq. (5.27) with Eq. (5.20), we see, except for the factor ZiZ^Z^^Z^^, the both identities are 
identical to each other. From this observation, we deduce 


Zi _ Zi 

%~Tz 

This is the S-T identity which coincides with the one given in the massless QCD [19, 20]. 


(5.28) 
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VI. GLUON FOUR-LINE VERTEX 


By the similar procedure as deriving Eqs. (5.1) and (5.3), the W-T identity obeyed by the gluon four-point Green 
function may be derived by differentiating Eq. (4.1) with respect to the sources J\{z) and Jf{u). The result 

represented in the operator form is as follows 


=< Q+\T*[C-{x)bl^{y)C^{y)Al{z)Ai{u)]\Q- > 

-b < 0+\T*[C^^{x)Al{y)bl^{z)C\z)Ai{u)]\G- > 

+ < Q+\T*[C^{x)Al{y)Al{z)bi%u)C%u)]\Q- > 

where 

G;itix,y,z,u) =< 0+|T[i;(x)i^(y)i^(z)i^(«)]|0- > (6.2) 

is the gluon four-point Green function. The accompanying ghost equation may be obtained by differentiating Eq. 
(4.2) with respect to the sources Jx{z) and J^iu). The result is 

< QbTbDlbx)C^{x)C\y)Al{z)Ai{u)]\Q- > .g 3^ 

+a^G‘^\‘l^{x,y,z,u) = -5‘^^5'^{x-y)Df^{z-u) ^ ' ’ 

where 

G^\i{x,y,z,u) =< 0+|r[G“(x)G^(y)i=^(z)i^(w)]|0- > (6.4) 


is the four-point gluon-ghost particle Green function. Differentiation of Eq. (6.1) with respect to the coordinates y, 
z and u and use of Eq. (6.3) lead to 


idl^d^d^dlG’lt^^.ix, y, z, u) = - y)d^dlD-^,{z - u) 

+6‘^‘=S'^{x - z)d'^dlDli{y - u) + - u)d'^d^Dl\{y - z) 

+<J^{d^,dlG'>‘^li{y, X, z, u) + d-dlG--ti{z, X, y, u) 
+d';dlG<^‘^^^l{u,x,y,z)}. 


It is noted that the four-point Green functions appearing in the above equations are unconnected. Their decompo¬ 
sitions to connected Green functions are not difficult to be found by making use of the relation between the generating 
functionals Z for the full Green functions and W for the connected Green functions as written in Eq. (3.9). The 
result is 

Gf^^tix, y, Z, u) = y, z, - Df,{x - y)Dli{z - u) 

-U-(x - z)Dli{y -u)- Df^ix - u)Dtliy - z) 

and 

G“*'^^(x, y, z, u) = G^^liix, y, z, u), - A^\x - y)Dt{z - u). 


( 6 . 6 ) 

(6.7) 


The first terms marked by the subscript ”c” in Eqs. ( 6.6) and (6.7) are connected Green functions. When inserting 
Eqs. (6.6) and (6.7) into Eq. (6.5) and using the W-T identity in Eq. (4.5), one may find 


(x, y, z, u), = x, y, z), 

+a);a;G““(z, x, y, u), + x, z, u) j. 


( 6 . 8 ) 


This is the W-T identity satisfied by the connected four-point Green functions. In the Landau gauge, we have 

d^dp^dlG^lt ix, y, U w)e = 0 (6.9) 


which shows the transversity of the Green function. 

The W-T identity for the four-line proper gluon vertex may be derived from Eq. (6.8) with the help of the following 
one-particle-irreducible decompositions of the connected Green functions which can easily be found by the standard 
procedure [13-16]. 
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G'IZ'ti^l,X2,X3,X4)c 

= I Ud‘^y^D1,p{xl - yi)Dl’>',{x2 - y2)T'^,b'c'd' (yi, 2/2, ys, ^ 4 ) 

x^A'Kys - X3)Df,‘^{y4 - X 4 ) 

+*/ n d‘^yid‘^Zi{Dii,{xi - yi)D^X’{^2. - y2)K-b-eiyi^y2,y5) 

i—1 

X^pp'(y3 - Zl)K'c'd' (^l> ^ 2 , Z3)^a'a(^2 - X 3 )Df,f{z 3 - X 4 ) 
+ ■ 0 ““'/ {xi - yi)Dl{, {X3 - y2)K'c'e (yi > y2, y3)^pp' (y3 - Zi) 

^ ^K’b’d' izi,Z 2 , Z3)D^Xiz2 - X2)DXr{z3 - X 4 ) 
+DXix-2 - yi)DXi.^3 - y2)rbv\'’(yi,y2,y3)£'pp'(y3 - 2 = 1 ) 
^K'a'd' {zi, Z 2 , Z 3 )DX^{z 2 “ Xi)Df 4 {z 3 - X 4 )} 


( 6 . 10 ) 


and 


GXri^l>^2,X3,X4)c 

= I n d'‘y*A““'(xi - i/i)r„,t,/c''j’''(yi,y2,y3,y4)A^'^(j/2 - X 2 ) 

^ XddX(y3-X3)DX(y4-X4) 

+if n d\*(i‘‘ 2 ;i{A““'(xi - yi)r^,gd'(yi,y2,y3)A“'(y2 - zi) 

xOXiyS - X4)^e'b'c'iZl,Z2, Z3)A^'^^{z2 - X2)DXiz3 “ 2 : 3 ) 
+A““'(xi - yi)r^,g^,'(j/i,y2,y3)A“ (y2 - zi)IlA''A(y3 - 2 ^ 3 ) 

X'^e'Jd'(Zl,Z 2 , Z3)A^'^{z 2 - X2)Df,f{z3 - X 4 ) 

+A““'(xi - yi)r^,b(’e(yi,y2,y3)A'''''(j/2 - X2)DZ,{y3 - zi) 


^X'c'd' {zi,Z 2 , Z3)DX{z2 - X 3 )DX{z 3 - X 4 )} 

where r“^‘ 3 ((.(xi,X 2 ,X 3 ,X 4 ) is the four-line gluon proper vertex and X 2 ,X 3 , 0 : 4 ) is the four-line ghost-gluon 


( 6 . 11 ) 


proper vertex. They are defined as [13-16] 


rpjAT(^l: 2^2,2;3, X4) — SA‘^i^(xi)SA0‘'(x2)SA‘^^(x3)SA<^^(x4) |2=0, 

r At ( 2^1 J 2 ^ 2 , 2 : 3 , 2 : 4 ) = iSC‘^(xi)SC’>(x2)SA<=>'(x3)SA'i^(x4)^'^=°' 


( 6 . 12 ) 


When substituting Eqs. (6.10) and (6.11) into Eq. ( 6 . 8 ) and transforming Eq. ( 6 . 8 ) into the momentum space, one 
can find the following identity satisfied by the four-line proper gluon vertex 


X^kXlKtiki,k2,k3,k4) = 1 :^ 


^3 
+'S 


fiuXr 

a c d b 
ki k3 k4 k 2 




a d b c 
ki k4 k 2 k3 


(6.13) 




where 

abed 

ki k 2 k3 k4 

= -ikX^^KXa{kuk2, -(fci + fe))i^:;(fci + k2)kXlXtX-^^3 + fc4), A:3, A:4) 

+ ^-^x{kl)x{kl)[ikXl^Xiik2, k4,k3, k4) (6.14) 

-A*'-(fe, A:i, -(A:i + fe))i?:;(fci + k2)k^klAfX{-{k3 + ^ 4 ), ^ 3 , k4) 

-klAX{k2,-{k2 + k4), k4)AG{k2 + k4)klA^"X-iki + k3), kiM) 

-kXXik2,-{k2 + k3), k3)AG{k2 + k3)klAf-i{-{ki + k4), ki, ^ 4 )]. 

The second and third terms in Eq .(6.13) can be written out from Eq. (6.14) through cyclic permutations. In the 
above, we have defined 


T;Xrikuk2,k3, k4) = {2 t:Y 6\Y.U h)AX^Xki, ^ 2 ,^ 3 , ^^ 4 ), 

TXi{kuk2,k3,k4) = {27TyS\Y:tih)^Xi{ki,k2,k3,k4). 


15 



In the lowest order approximation, we have checked that except for the first term in Eq. (6.14) which was encountered 
in the massless theory, the remaining mass-dependent terms are cancelled out with the corresponding terms contained 
in the second and third terms in Eq. (6.13). Therefore, the identity in Eq. (6.13) leads to a result in the lowest order 
approximation which is consistent with the Feynman rule. 

The renormalization of the four-line vertices is similar to that for the three-line vertices. From the definitions given 
in Eqs. (6.12), (6.15) and (5.18), it is clearly seen that the four-line vertices should be renormalized in such a manner 

h, h) = h, h), 

fc 2 , ks, ki) = zYzY^YYAkiA2,k3, h). ^ ^ 

On inserting these relations into Eqs. (6.13) and (6.14), one can obtain a renormalized identity similar to Eq. (5.20), 
that is 






k^k!^kYlAYYirikl, k2, kM = '^R ( I 

a c d b 
ki fca ki k 2 


1 '>-2 
a d b 
ki ki k 2 



where 


f a b c d \ 

V ki k2 k, ki ) 

= -ik^ik^YtlAkiM, -(fci + k2))DAYki + k2)k^klAkA;^^AYk3 + kAMAi) 

+ ^Xit(*i)Xit(fci)[*fc3^4Afl'='l.,-(fc2, ki, ks, ki) 

-AYAk2Ai,-iki + k2))DAA{ki + k2)k^klAllt{-{k^ + ki), k,,ki) 
-klAYrik2,-ik2 + ki), ki)A^jl(k2 + ki)k^AYli-iki + k^), ki, k^ 
-k^AYx{k2,-{k2 + kA,ki)Ai (*2 + kAklY^A-iki + ki), ki,ki)]. 


(6.17) 


(6.18) 


We can also define vertices AYxriki, k2, ks, ki) and A'^Yriki, k2, ks, ki) from the vertices A^Aiki, k2, k^, ki) and 
A°'Y{ki,k 2 ,k 3 ,ki) by taking out the coupling constant squared, respectively. The renormalization of these vertices 
are usually defined by [13-17] 


Y^YAkiA 2 , ks, ki) = Z 4 ^~AiAYAkiA 2 , ks, ki), 

AYAki,k2,k3, ki) = Z^Aj^Yriki, k 2 , k^, ki). ^ ’ 

where Zi and Zi are the renormalization constants of the four-line gluon and ghost-gluon vertices respectively. 
Obviously, the identity in Eqs. (6.13) and (6.14) remains formally unchanged if we replace all the vertices in the 
identity with the ones . Substituting Eqs. (5.24), (6.19), (4.26) and (4.27) into such an identity, one may write a 
renormalized identity similar to Eq. (5.25), that is 


LA* I,A l,T A abed 
^1 ^2 ^3 ^ r 


'■3 


(fci, k 2 ,k 3 , ki) = 'I’r 


ki 




a c d b 
ki fcs ki k 2 


+ 4*^ 


a d 
ki ki 


b c d \ 

k 2 ks ki J 



( 6 . 20 ) 


where 


_ Z4Z3 


~ f a b c d \ 

^^\krk2k3ki) 

{-tk'^k!^AYAAkiA2,-iki + k2))DA!:f{ki + k2)k^klAl;li-{k3 + ki), k^, ki 
+ ^-^XR{kl)xR{kl){^jkYYYY{k2M,k3, ki) 

|^A^“ (fc2, -(fci + k2))DAYki + k2)klklAl;l{-ik3 + ki),k3, ki) 
-^[klAYAk 2 ,-{k 2 + ki), ki)A)f {k 2 + ki)kfAl\{-{ki + ks), ki,k3) 
+k^AYxik2,-{k2 + ki), k3)A^jk (fe + A:3)fcJA/“^(-(fci + ^^ 4 ), A:i, ^ 4 )]}. 


)} 


( 6 . 21 ) 
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Multiplying the both sides of Eqs. (6.20) and (6.21) by according to the relations given in Eqs. (5.26) and in the 
following 


k2, fcs, ki), 

fe) ^3, k4) = gRJ^R°’^Xr{kl, ^2, A:3, ^ 4 ), 


( 6 . 22 ) 


we have an identity which is of the same form as the identity in Eqs. (6.20) and (6.21) except that the vertices in 
Eqs. (6.20) and (6.21) are all replaced by the vertices K^r. Comparing this identity with that written in Eqs. (6.17) 
and (6.18), one may find 


which lead to 


Z 3 Z 4 _ 

’ Z 3 Z 4 


y Z/^Zz _ ^ 4^3 

’ ZiZi “ ’ ZzZl 


1 


Zi Z^ Z\ Z'^ z^ 

Zz~ Zz~ Z4~ Zz~ Z 4 


(6.23) 


(6.24) 


This just is the S-T identity which is consistent with that given in Refs. (19) and (20) for the massless QCD. 


VII. QUARK-GLUON VERTEX AND QUARK PROPAGATOR 


This section is used to derive the W-T identity for quark-gluon vertex and discuss its renormalization. First we 
derive a W-T identity satisfied by the quark-antiquark-gluon three-point Green function. This identity can easily be 
derived by differentiating the W-T identity in Eq. (3.8) or (3.10) with respect to the sources ri{y) and ^(x) and 

then setting all the sources to be zero. The result written in the operator form is as follows 



d^G^ix, y, z) = iaglG^ix, y, y, z)T^ - T'’G^“(x, y, x, z)] 

(7.1) 

where 

Gl{x,y,z) = ^ 

0 + 1 ^{x)^{y)A{z) 1 0 ^ 

(7.2) 

is the quark-gluon 

three-point Green function, 




GJ“(x,2/,z) = (0+ 

1 ${x)^{y)C^{y)G {z) 1 0’^ 

(7.3) 

and 

G 2 °'{xry,z) = ( 0 + 

1 i^{x)^{y)C\x)d {z) 1 0~\ 

(7.4) 


are the quark-ghost particle mixed Green functions. The Green functions in Eqs. (7.3) and (7.4) are connected 
because a quark field and a ghost field are of a common coordinate. 

The W-T identity for quark-gluon vertex can be derived from Eq. (7.1) with the help of one-particle irreducible 
decompositions of the Green functions shown in Eqs. (7.2)-(7.4). The decompositions can easily be obtained by the 
standard procedure [13-16]. The results are given in the following. 

y,z) = J d’^x'd’^y'd'^z'iSpix - x')T^''{x', y', z')iSF{y' - y)iDl°'^{z' - z) (7.5) 

where D’l°-^{z' — z) is the gluon propagator defined in Eq. (4.6), 


is the quark propagator and 


iSpix — x') = (0'*' I 4’i.x)4){x') I 0 


r OL' / / / /\ 

(x ,y , z ) = —=— 

^ ^ ' iXnl,Or./ 




i6il}{x')5'4i{y')5A^^{z') 


|j=o 


(7.6) 

(7.7) 
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is the quark-gluon proper vertex. 


G\°'{x,y,z) = j d'^x'd!^z'S f{x — ,y, — z) 

where A“(z' — z) is the ghost particle propagator dehned in Eq. (4.7) and 

y, z') = J d'^ud'^vA^'^iy - u)T‘'‘^{x', v, u, z')Sf{v - y) 


in which 


V, u, z') = i 


6^r 


6ip{x')6'ijj{v)6C {u)6C‘^{z' 


|j=o 


is the quark-ghost vertex. Similarly, 



d^y'd'^z'Y^ix, y', z')SF{y' 






where 


y\ z') = J d-^ud^^ySFix - u)A^‘^{x - v)T‘^%u, y', 


(7.8) 

(7.9) 

(7.10) 

(7.11) 

(7.12) 


On substituting Eqs. (7.5), (7.8) and (7.11) into Eq. (7.1) and then transform Eq. (7.1) into the momentum space, 
we have 


SFip)T^''iP,q,k)SF{q)k^^D-iik) 

= -iag[SF{p)lKp, q, k) - jUp, 9, k)SF{q)]A°-^{k)] 


(7.13) 


where we have defined 

jUP,q,k)=7fip,q,k)T\ 

j^{p,q,k)=T\l-{p,q,k). 


(7.14) 


Considering that the vertex functions r^''{p,q,k) and 'yf{p,q,k) {i = 1,2) contain a common delta-function repre¬ 
senting the energy-momentum conservation, we may set 

r“^(p, q, k) = {2Tr)'^S'^{p -q + fc)A“'"(p, q, k), , . 

7f{p,q,k) = {2Tr)‘^S'^{p-q + k)-ff{p,q,k) 


where h°'^^{p, q, k) and jf{p, q, k) are the new vertex functions in which k = q — p. Noticing the above relations and 
the expressions of gluon and ghost particle propagators as given in Eqs. (4.11) and (4.22), the W-T identity in Eq. 
(7.13) can be rewritten via the functions A°'^{p, q, k) and 7 “(p, g, k) in the form 

k^A°-^^{p,q,k) = igx{k‘^)[SF^{p)72{P^q,k) -JiiP,q,k)Sp'^{q)] (7.16) 

where x(fc^) '^as defined in Eq. (5.13). 

Let us turn to discuss the renormalized form of the above W-T identity. It is well-known that the quark propagator 
can be expressed in the form 


SfIp) 


1 

p—m — E(p) -I- ie 


(7.17) 


where p = j^Pp and E(p) denotes the quark self-energy. The above expression can easily be derived from the Dyson 
equation [29]. Usually, the quark propagator is renormalized in such a fashion 

SFip) = Z^S^ip) (7.18) 


which implies 


V’(x) = ^/^ipRix), -tpix) = ^/^■^l)R{x). 


(7.19) 
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From the relations in Eqs. (7.19) and (5.18), it is clearly seen that the vertex defined in Eq. (7.7) is renormalized as 

r“^(x,y,z) = Z-^ZpTl^{x,y,z) (7.20) 

which leads to 

A“^(p, q, k) = Z-^zpAl^{p, q, k) (7.21) 

The functions 7 “ (p, q, k) are, in general, divergent in the perturbative calculation. These functions are assumed to be 
renormalized in such a manner 


7 ?(p, q, k) = Z^ q, k). (7.22) 

where Z^ is the renormalization constant of the functions jf{p,q,k). Based on the relations in Eqs. (5.22), (7.18), 
(7.21), and (7.22), Eq. (7.16) can be represented in terms of the renormalized quantities 

k^Lki‘^R{p,q,k) = i9RXR{k'^)[SF~^{p)l2RiP,q,k)-JiRip,q,k)S^~'^{q)] (7.23) 

where is the renormalized coupling constant defined by 

gR = %zlz-^g (7.24) 

It is well-known that 

gR = %zlz-^g (7.25) 

where Zi is the ghost vertex renormalization constant as defined in Eq. (5.24). The relation in Eq. (7.25) ordinarily 
is determined from the renormalization of S-matrix elements. In comparison of Eq. (7.24) with Eq. (7.25), we see 

Z^ = Zi (7.26) 

which means that the functions 7 “ (p, q, k) are renormalized in the same way as for the ghost vertex. 

In the conventional discussion of the vertex renormalization, one considers such a vertex denoted by h°'^{p,q,k) 
that it is defined from A°'^^{p,q,k) by taking out a coupling constant. Obviously, the W-T identity obeyed by the 
A°'>^{p, q, k) can be written out from (7.16) by taking away the coupling constant on the RHS of Eq. (7.16), that is 

k^A^^{p,q,k) = ix{k'^)[Sp^{p)x2{P,q,k) -xl{p,q-,k)Sp^{q)]. (7.27) 

The renormalization of the vertex A°‘^{p,q, k) usually is defined by 

A<^^^{p,q,k) = Zp^A%^{p,q,k) (7.28) 

where Zp is the quark-gluon vertex renormalization constant. When Eqs. (5.22), (7.18), (7.22) and (7.28) are inserted 
into Eq. (7.27) and then multiplying the both sides of Eq. (7.27) with a renormalized coupling constant, we arrive at 

k^,A°,^{p,q,k) = iZpZ2,Z:^^Z-^gRXR{k'^)[Sp~^{p)l2R{P,q,k) -ilR{p,q,k)S^~^{q)\ (7.29) 

where 

J^°R{p,q-,k) = gRA°,^{p,q,k). (7.30) 

In comparison of Eq. (7.29) with Eq. (7.23) and considering the equality in Eq. (7.26), we find, the following identity 
must hold 


Zp _ Zi 

Combining the relations in Eqs. (5.28), (6.24) and (7.31), we have 

Zp Z\ Z\ Z/^ 
'Z‘2 Z-i Zi 


(7.31) 


(7.32) 


This just is the well-known S-T identity. This identity was obtained from the massless QCD and now, as has just 
been proved, it also holds for the massive QCD. 
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VIII. EFFECTIVE COUPLING CONSTANT AND GLUON MASS 


This section and the next section are used to perform one-loop renormalization of the massive QCD by using the 
renormalization group approach. As argued in our previous paper [30-32], when the renormalization is carried out in 
the mass-dependent momentum space subtraction scheme, the solutions to the RGBs satished by renormalized wave 
functions, propagators and vertices can be uniquely determined by the boundary conditions of the renormalized wave 
functions, propagators and vertices. In this case, an exact S-matrix element can be written in the form as given in 
the tree-diagram approximation provided that the coupling constant and particle masses in the matrix element are 
replaced by their effective (running) ones which are given by solving their renormalization group equations. Therefore, 
the task of renormalization is reduced to hnd the solutions of the RGBs for the renormalized coupling constant and 
particle masses. Suppose Br is a renormalized quantity. In the multiplicative renormalization, it is related to the 
unrenormalized one F in such a way 


F = ZfFr (8.1) 

where Zp is the renormalization constant of F. The Zp and Br are all functions of the renormalization point /r = /ioe* 
where hq is a fixed renormalization point corresponding the zero value of the group parameter t. Differentiating Bq. 
(8.1) with respect to /i and noticing that the F is independent of /i, we immediately obtain a renormalization group 
equation (RGB) satisfied by the function Fr [21-23] 

dF 

+ IfFr = 0 ( 8 . 2 ) 

dll 

where 7 f is the anomalous dimension defined by 


7f = 11 -^ InZp- (8.3) 

dll 

Since the renormalization constant is dimensionless, the anomalous dimension can only depend on the ratio fd = 
where tur denotes a, renormalized mass and 7 r= " 1 p{gR,ld) in which qr is the renormalized coupling constant and 
depends on /i. Since the renormalization point is a momentum taken to subtract the divergence, we may set g, = /rgA 
where A = e* which will be taken to be the same as in the scaling transformation of momentum p = poX. In the 
above, po and po are the fixed renormalization point and momentum respectively. When we set F to be the coupling 
constant g and noticing = A^, one can write from Bq. (8.2) the RGB for the renormalized coupling constant 


A 


dgniX) 

dX 


+ lg{X)gR{X) = 0 


(8.4) 


with 


7. = P-1„2,. 


According to the definition in Bq. (8.1) and the relation in Bq. (7.25), we may take. 


(8.5) 


Zg = (8.6) 

to calculate the anomalous dimension. As denoted in Bqs. (4.25) and (5.24), the renormalization constants Z 3 , Z 3 and 
Zi are determined by the gluon self-energy, the ghost article self-energy and the ghost vertex correction, respectively. 
At one-loop level, the gluon self-energy is depicted in Figs, (la)-(ld), the ghost article self-energy is shown in Fig. 
(2) and the ghost vertex correction is represented in Figs. (3a) and (3b). According to the Feynman rules which are 
the same as those for the massless QCD [16] except that the gluon propagator and the ghost particle one are now 
given in Bqs. (4.14) and (4.10), the expressions of the self-energies and the vertex correction are easily written out. 
For the gluon one-loop self-energy denoted by —zn“]((A:), one can write 

4 

= (8.7) 
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where (k) represent the self-energies given in turn by Figs.(la)-(ld). They are separately represented 

in the following: 


n 


(l)a 6 

liv 


{k) — f + 2A:)p g}.p(2l + k)^ 

~k9p^^Q‘ ~ k)x\[9i'p'Q- ~ k)y — g\ip'{2l + k)i, + gy^il + 2k)pi], 


( 8 . 8 ) 




d^l {I k)^li, 

(27r)4 [{I + kY — + ie] [l'^ — M'^ + is] ’ 


(8.9) 


and 


(k) — id 5g J ^^2 _ jy^2 _|_ jg) i9pv9xp 9pp9xv) 


-p|-(4)a6 


(fc) 


_ „-xo61„2 f d-^l 1 

- -10 2ff J JW)^ [i^i^ky^^m‘^+ie][l‘^-m‘^+ie] 
xTr[ 7 p(l - k -I- to)7i,(1 -I- m)] 


( 8 . 10 ) 


( 8 . 11 ) 


where 1 k =j^k\. In the above, J“cdj6cd _ ^gab Tr{T°'TY = have been considered. It should 

be noted that in writing Eqs. (8.8)-(8.10), we choose to work in the Feynman gauge for simplicity. This choice is 
based on the fact that the massive QCD has been proved to be an unitary theory [18], that is to say, the S-matrix 
elements evaluated from the massive QCD are independent of gauge parameter. Therefore, we are allowed to choose 
a convenient gauge in the calculation. From Eqs. (8.8)-(8.11), it is clearly seen that 


UfYk) = S^'^UpYk) = ^“'^n«(A:). (8.12) 

i=l 


By the dimensional regularization approach [33-37], the divergent integrals over I in Eqs. (8.8)-(8.11) can be regularized 
in a n-dimensional space and easily calculated. The results are 


Ilpj(fc) — 2 ( 4^)2 /o '^^£[fc2x(a:-l) + M2]s { 5 /ii/[lla;(x 1 ) 

+5)fc^ -I- 9M^] -I- 2[5a:(x — 1) — l]kpk^}, 


(8.13) 


n|fj(fc) = fg da^ ,[fc.x(x-\)+M 2 )e {[fe"^(a^ - 1) 

+M^]gp^ + 2x(x - l)kpk!,}, 


(8.14) 


and 



9g2 m 2 



4g^ Px{x-1) kpki, 

(47r)2 _/g ^£[kP‘x(x — 1) +tv?Y kP' 


(8.15) 


(8.16) 


where e = 2 — ^ > 0 when n ^ 4. In Eqs. (8.13)-(8.16), except for the e in the factor I/e[fc2a;(a; — 1) -I- M^]® and 

l/£[k‘^x{x — I) -f rnPY, we have set e ^ 0 in the other factors and terms by the consideration that this operation does 
not affect the calculated result of the anomalous dimension. According to the decomposition shown in Eqs. (4.15) 
and (4.16) and noticing g^^, = T’Y' + V'Y', it is easy to get the transverse part of Wp^^k) from Eqs. (8.13)-(8.16) and 
furthermore, based on the decomposition denoted in Eq. (4.20), the functions ni(A:2) and 112(^2) can be written out. 
The results are 


ni(fc 2 ) 


r 15[2a:(a:-1)-kl] ^ 4a:(x - 1) 

{AttY Jq ^ 2£:[A:2a;(x — 1)-I-M2]e £[fc2a:(x — 1 )-I-to2]® 


(8.17) 


and 
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(47r)2^/ '^'^e[k‘^x{x-l)+M^Y 2£M2e^- 

It is clear that the both functions ni(fc2) and Ii 2 (k^) are divergent in the four-dimensional space-time. When the 
divergences are subtracted in the mass-dependent momentum space subtraction scheme [34-37], in accordance with 
the definition in Eq. (4.25), we immediately obtain from the expression in Eq. (8.17) the one-loop renormalization 
constant as follows 


Z3 = l-ni(/x 2 ) 

15[2a;(a:-l) + l] ^ 4x{x-l) ^ (8.19) 

( 47 r )2 Jo “•^l 2 £[/i 2 a;(a;_i)+M 2 ]s s[n‘^x(x-l)+mj]^ S' 


Next, we turn to the ghost particle one-loop self-energy denoted by From Fig. (2), in Feynman gauge, 

one can write 




dH 


q-{q-l) 


(27r)‘* [{q — ly — M2 -|- ie\[P — M'^ + ie] 


By the dimensional regularization, it is easy to get 

n'^\q) = 5'^\^h{q^) 

where 


(l{q^) = 


r 


dx 


3(a: — 1) 


(47r)2 Jp e[g2a;(a; - 1)-I-M2]e ■ 


( 8 . 20 ) 


( 8 . 21 ) 


( 8 . 22 ) 


According to the definition given in Eq. (4.25) and the above expression , the one-loop renormalization constant of 
ghost particle propagator is of the form 


^3 = 1 - niixY = 1 - 




dx 


3(a;- 1) 


( 47 r )2 7 p e[fj,'^x{x — 1 ) + M^Y 


(8.23) 


Now, let us discuss the ghost vertex renormalization. In the one-loop approximation, the vertex defined by 
extracting out a coupling constant is expressed as 


Ar(p, q) = r^px + ATxip, q) + (p, q) 


(8.24) 


where the first term is the bare vertex, the second and the third terms stand for the one-loop vertex corrections shown 
in Figs. (3a) and (3b) respectively. In the Feynman gauge, the vertex corrections are expressed as 


AfYip^q) = 


dH 


p ■ iq-l)ip-l)\ 


{2ttY [P — M2 -I- ie] [(p — ly — M2 -|- ie] [(g — ly — M2 -|- ie] 


and 


A abc f ^ \ • nabc^ ^2 

2A (P: q) = */ 2^ 


dH 


I - [p- q - l)px - pHqx+p - {2q-p + l)lx 


(27r)4 [P — M2 -I- ie] [(p — ly — M2 -|- ie] [(g — i)2 — M2 -|- ie] 


(8.25) 


(8.26) 


where f°-<^d.j;ebfj;dfc _ _^jabc jpggg noted. By employing the dimensional regularization to compute the above 
integrals, it is not difficult to get 


Afy{j),q) = [ dx f dy{j^ - -^[pxAi{p,q) + qxBi{p,q)] - ^pa} (8.27) 

2^ fWj Jo Jo '^xy o 

where 

Qxy = p^xy{xy -1) + q'^lix - lYy + {x - l)]y - 2p • qx{x - l)p2 -h M'^, 

Ai{p,q) = {p-q[l + {x- l)y\ -p^xy}{l - xy)y, (8.28) 

■Bi(p, q) = {p- g[l + {x- l)p] - p‘^xy}{x - l)p2 
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and 




2 (47r)2 


/o 




fj/PA 


1 3 

+ t;—[pa^ 2 (p, g) + q\B2{p, g)] - -p\} 

^XV o 


(8.29) 


where 


M{p,q) = {p^{2xy - x'^y'^ - 1 ) - q‘^[{x - l)y - l](a;- l)y 
+P ■ g[2 - (3a: - 2)y + 2x{x - l)y'^]}y, 

B 2 {p,q) = [p ■ q{x - 1 ) -p^x\y'^. 


(8.30) 


The divergences in the both vertices K‘^^{p,q) and K‘^^{j},q) may be subtracted at the renormalization point p^ = 
q^ = p? which implies k = p — q = 0, being consistent with the momentum conservation held at the vertices. Upon 
substituting Eqs. (8.27) and (8.29) in Eq. (8.24), at the renormalization point, one can get 


Af=(p,g) + Li) = Z^^r^p^ 


nabc 


— 1 xabc^ 


(8.31) 


where 


Zi = 1 - Li = 1 - 


3g^ 

(47r)2 7o ^'^^e[p^x{x-l) + M^Y 


dx{- 


x'^{x — l)pY 1 

p^x{x — 1) + M2 4 


(8.32) 


which is the one-loop renormalization constant of the ghost vertex. 

Now we are ready to calculate the anomalous dimension 7g(A). Substituting the expressions in Eqs. ( 8 . 6 ), (8.19), 
(8.23) and (8.32) into Eq. (8.5), it is easy to find an analytical expression of the anomalous dimension 7g(A). When 
we set ^ j and ^ ^ with defining /? = ^ and ^ (here we have set po = A), the expression of 7g(A), in 

the approximation of order g^, is given by 


9r 


7,(A) = In Z. - In Z 3 - In Z,] = 


F(A) 


(8.33) 


where 


IZCAA _ 19 lap I SA"' /'o lop* A" A ap" 

IM - — - ^2- + 2(Aa-4P^) - - ^2- - 

X In ■ 




-iE[i 


^Pi 


VlPi 


: In 


X-^X^-Apl 


(8.34) 


A+^A2-4/32 "i ^ A3^A2-4p2 A-|-^A2-4p2 


in which Nf denotes the number of quark flavors. We would like to note that the fixed renormalization point A in /3 
and Pi can be taken arbitrarily. For example, the A may be chosen to be the mass of the quark of ASj -th flavor. In 
this case, /? = M/iriNf and pi = mi/ruNf In practice, the A will be treated as a scaling parameter of renormalization. 

With the 7g(A) given above, the equation in Eq. ( 8 . 4 ) can be solved to give the effective coupling constant as 
follows 


afl(A) = 


CtR 


l + ffG(A) 


where Q;fi(A) = g^(A)/ 47 r, aR = ^/{(l) and 


(8.35) 


/•A j\ 1 

G(A) = -F{X) = pi{X)-Ml)--YYpl{X)-pl{l)] 


(8.36) 


in which 


10/32 Ma 2 -4/32 
pi{X) = [(19-^)- 


A 2 


4 A 


3 A 1 A + VA 2 - 4/32 

4 a/A 2 -4/32^ A - a/A 2 - 4/32 


5/32 


(8.37) 


.>.(A)-(l + ^)-j-■" A - yA» - 4pi ~ V- 


(8.38) 
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(8.39) 


and v?i(l) = ViW |a=i, ^2^ — ^2.{^) L=i- large momentum limit (A ^ oo), we have 

G(A) = (ll-^iV/)lnA. 

This just is the result for massless QCD which was obtained previously in the minimal subtraction scheme [38-40]. 
It should be noted that the expressions in Eqs. (8.34), (8.37) and (8.38) are obtained at the timelike subtraction 
point where the A is a real variable. We may also take spacelike momentum subtraction. For this kind of subtraction, 
corresponding to /r —> ifi, the variable A in Eqs. (8.34), (8.37) and (8.38) should be replaced by iX where A is still a 
real variable. It is easy to see that the function in Eq. (8.39) is the same for the both subtractions. 

The behavior of the function aR{X) is graphically described in Figs. (4)-(6). Figs. (4) and (5) represent respectively 
the effective coupling constants obtained at the timelike subtraction point and the spacelike subtraction point, where 
we take the flavor Nf = 3 as an illustration. For comparison, we show in each of the figures three effective coupling 
constants which are obtained by the massive QCD, the massless QCD and the minimal subtraction, respectively. 
These effective coupling constants are respectively represented by the solid, dashed and dotted lines in the figures. 
To exhibit the dependence of the effective coupling constant on the number of quark flavors, in Fig. (6), we show 
three effective coupling constants given by the timelike momentum subtraction. The solid, dashed and dotted lines 
in the figure represent the effective coupling constants obtained by taking Nf = 2, 3 and 4, respectively. In our test, 
we find that the behavior of the Q;ij(A) sensitively depends on the choice of the constant an, the scaling parameter A, 
gluon mass M/j and the quark masses m/j. Certainly, the parameters a^, A, and m/j should be determined by 
fitting to experimental data. In our calculation, as an illustration, we take the quark masses to be constituent quark 
ones. The masses of up, down, strange and charm quarks are taken to be = nid = 350MeV, nig = 500MeV and 
nic = ISOOMeF. For the other parameters, we take = 0.2, = QOOMeV, A = 500MeV in Figs. (4) and (5) and 

A = 1500MeV in Fig. (6). From Fig. (4) it is seen that the effective coupling constant given by the massive QCD 
is an analytical function with a maximum at A = 1.346, the effective coupling constant given by the massless QCD 
is also an analytical function with a peak around A = 1, but the effective coupling constant given by the minimal 
subtraction has a singularity at A = 0.1746. Fig. ( 5 ) indicates that the effective coupling constant given by the 
massive QCD (where gluon mass Mr = 600Mefo), similar to the one given by the minimal subtraction, has a Landau 
pole at A = 0.1845 which implies that the coupling constant is not applicable in the region A < 0.1845. However, 
if the gluon mass is taken to be Mr < 425.75MeF, we find, the Landau pole disappears and the effective coupling 
constant, analogous to the effective coupling constant given by the massless QCD, becomes a smooth function in the 
whole region of momentum as illustrated by the dotted-dashed line in Fig. (5) which represents the effective coupling 
constant given by taking Mu = 425.75MeF. From Fig. (6) it is clear to see that in the low and intermediate region 
of momentum, the larger the flavor number Nf^ the smaller the maximum of the effective coupling constant is and 
the positions of the maxima for different coupling constants are different from one another. This property of the 
effective coupling constants would give a notable effect on the theoretical hadron spectrum because for the calculation 
of hadron spectrum, the quarks in a hadron are assumed to move not too fast as suggested in the nonperturbative 
quark potential model, therefore, the behavior of the effective coupling constant in the low and intermediate domain 
of energy would play a dominate role In this case. 

Let us proceed to derive the one-loop effective gluon mass. Setting Fr = Mr in Eq. (8.2), we have the RGE for 
the renormalized gluon mass 


A 


dMR{X) 

dX 


+ 1m{X)Mr{X) = 0 


(8.40) 


where 


7m(A) = IuZm- 
dfx 


(8.41) 


From the last equality in Eq. (4.25) and Eqs. (8.17) and (8.18), in the approximation of order we can write 


Zm — 1 + 5[ni(M^) + 112 (/^^)] 

_ , _ r j r 3[10a:(a:-l)-|-13] , 2x(x-l) i _ 27 1 (8.42) 

(47r)2 ^ Jo ^4€[k‘^x{x — l)-\-M‘^]^ e[k'^x(x — l)-\-m‘'^]^ ^ 

i—1 ^ 

On inserting Eq. (8.42) into Eq. (8.41) and completing the differentiation with respect to /r and the integration over 
cc, we find 
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(8.43) 


7M(Aj - 

Nf 


3/3^(13A^-10/3^) A--^A^-4/3' 


-|E[i 


i=l 


6p- 


12pf 


In 


-4/3' 
A-^A^IIip 


AV^"-4p? A+^A2-4p: 


A+V^='-4/3= 
]}• 


With this anomalous dimension, the RGE in Eq. (8.40) can be solved to give an effective gluon mass such that 

M^j(A) = (8.44) 

where Mu = Mfl(l) and 


SgW=J^ y7M(A). 


(8.45) 


In general, the coupling constant gji in Eq. (8.43) may be taken to be the effective one. If the coupling constant is 
taken to be the constant gji, the function Sg{X) can be explicitly represented as 


(A) = ^{¥^3(A) - ^3(1) - ^ E[^2(A) - V^*2(l)]} 


Nf 


i=l 


where 


■.3(A) = (17 - —1” + ir 

and (/ 32 (A) was given in Eq. (8.38). In the large momentum limit, 

^3(A) = ^(17-^iV/)lnA. 

Therefore, we have 


(8.46) 


(8.47) 


(8.48) 


lim Mfl(A) = 0 


(8.49) 


which exhibits the asymptotically free behavior. 

The behavior of the effective gluon mass M/j(A) may be discussed in the way similar to the discussion of the effective 
coupling constant. Here we only limit ourself to show in Fig. (7) the behavior of the effective gluon mass written 
in Eqs. (8.44)-(9.46) with taking Nf = 3 and the coupling constant being a constant. In Fig. (7), the solid line 
and the dashed line represent the effective gluon masses given by the timelike momentum subtraction and spacelike 
momentum subtraction, respectively. The solid line exhibits that for the timelike momentum, the effective gluon mass 
is an analytical function with a maximum at A = 1.233. When A goes to infinity, the M/{(A) tends to zero rather 
rapidly, while, when A goes to zero, the Mr{X) abruptly falls to zero. The dashed line tells us that for the spacelike 
momentum, the Mr{X) keeps a constant Mr in the region [0,1] of A, while when A ^ 00 , the Mr{X) smoothly tends 
to zero. 


IX. EFFECTIVE QUARK MASS 

Before deriving the one-loop effective quark mass, we need first to discuss the subtraction of the quark one-loop 
self-energy on the basis of the W-T identity represented in Eq. (7.16). For later convenience, the identity in Eq. 
(7.16) will be given in another form. Introducing new vertex functions q) and 7 “(p, q) defined by 

r“^(p, q, k) = (27r)4(54(p -q + k)igM^^{p, q) .q . 3 

lt{p,q,k) =-{2 t:Y5'^{p- q + k)^f{p,q) ^ 

where i = 1, 2 and 7 “ (p, q) = — 7 “(p, q) and considering k = q — p, Eq. (7.16) can be rewritten as 

(p-9)pA“^(p,g) =x(fc^)[5'^^(p)72(p,g)-7i(p,g)*S'^\g)]- (9.2) 
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From the perturbative calculation, it can be found that In the lowest order, we have 


q) = q) = 'r“. 


(9.3) 


In the one-loop approximation of order the quark-gluon vertex denoted by g) is contributed from the 

two diagrams in Figs. (8a) and (8b) whose expressions can easily be written out. The quark-ghost vertex functions 
7i^^°(Pi'Z) (* = 1;2) are contributed from Figs. (9a) and (9b) and can be represented as 

#“(p,9)=™,(p,g) (9.4) 


where 

Kiip,q) =i^g‘^ J (9.5) 

and 

K2{p,q) =i^g^ J -^^SF{l)l^D^Yq-^){p-lTHp-^)- (9-6) 

It is clear that the above functions are logarithmically divergent. In the one-loop approximation, the function xO^) 
can be written as x(fc^) = 1 ~ where the ghost particle one-loop self-energy was represented in Eq. 

(8.22). Thus, up to the order of g"^, with setting q) = T“A^^^(p, g), we can write 

A“(p,g)=T“[ 7 ^ + AW(p,g)] (9.7) 

and 

x(fc2)7“(p,g)=r“[l + /,(p,g)] (9.8) 

where 

h{p,q) = KYp.,q) (9.9) 


Upon substituting Eqs. (9.7) and (9.8) and the inverse of the quark propagator denoted in Eq. (7.17) into Eq. (9.2), 
then differentiating the both sides of Eq. (9.2) with respect to and finally setting g = p, in the order of we get 


K{P,P) 


d^jp) 

dp>^ 


(9.10) 


where 


= A^^’iPyP) - 1,J2{P,P) - (P - »n) 

+2^ l,=p(p-m) 


dl2ip,q) 


(9.11) 


here m is the i-th quark mass (hereafter the subscript of rrii is suppressed for simplicity). It is emphasized that at 
one-loop level, the both sides of Eq. (9.11) are of the order of g^. In the derivation of Eq. (9.11), the terms of orders 
higher than g^ have been neglected. The identity in Eq. (9.10) formally is the same as we met in QED. By the 
subtraction at p= g, the vertex K^(p,p) will be expressed in the form 

^piP,p) = Ljf, + A^ip) (9.12) 

where L is a divergent constant defined by 

L = A^(p,p) |p=^ (9.13) 

and A^ (p) is the finite part of A^ (p, p) satisfying the boundary condition 

A;(p) Ip=^= 0. (9.14) 
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On integrating the identity in Eq. (9.10) over the momentum and considering the expression in Eq. (9.12), we 
obtain 


S(p) = A+{p>- d)[B- C{p^)] 

(9.15) 

where 


A = ^id), 

(9.16) 

B = -L 

(9.17) 

and C{p^) is defined by 


f dp>^Al{p) = (p - p)C{p^). 
dpa 

(9.18) 

Clearly, the expression in Eq. (9.15) gives the subtraction version for the quark self-energy which is required by the 
W-T identity and correct at least in the approximation of order g^. With this subtraction, the quark propagator in 
Eq. (7.17) will be renormalized as 

Sf{p) = , . 

p-mR- Zr[p) 

(9.19) 

where Z 2 is the renormalization constant defined by 


= l-B, 

(9.20) 

ruR is the renormalized quark mass defined as 


vfiR = Z~^m 

(9.21) 

in which 


= 1-1- Z2[Am~^ + (1 - P'm~^)B], 

(9.22) 

Zm is the quark mass renormalization constant and J1r{p) is the finite correction 
boundary condition S/{(p)|p 2^^2 = 0. 

Now we are in a position to discuss the one-loop renormalization of quark mass, 
quark mass can be written from Eq. (8.2) by setting F = m, 

of the self-energy satisfying the 

The RGE for the renormalized 

A ^ ^ +7m(A)TOij(A) = 0 
dX 

(9.23) 

where 


7m(A) = \nZm- 

dfx 

(9.24) 


It is clear that to determine the one-loop renormalization constant Zm, we first need to determine the divergent 
constants A and B from the self-energy represented in Eq. (9.15). The one-loop self-energy denoted by —i'Z{p) can 
be written out from Fig. (10). In the Feynman gauge, it is 


VM = -■- 2 /" _ 7^(k-f-p + TO)7^ _ 

*3^7 (2^)4 [(fc + p)2 - m2 + is] (P - M2 + ie) 


(9.25) 


where k and p By making use of the dimensional regularization to calculate the above integral, it is 

found that 


t - l)p+2m _ 

3 (4^)2 7o e[p^x{x - 1) + w?x + M^{1 - a;)]" ' 


(9.26) 
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According to Eq. (9.16), we have 


A = Sfni I - = [' dx _ ~ l)M+2m _^ 

^ ^ 3 (47r)2 7g e[^'^x{x — 1) + w?x + M'^{1 — x)Y 

With the aid of the following formula 


(9.27) 


6® Jq [aa; + 5(1 — a:)]i+® ’ 


one can get from Eqs. (9.15) and (9.27) that 


(9.28) 


8 g- 
3 (47r)^ 


fo 


B = [E(p) - A](p-to) 

l)+m^fc+M^(l —ai)]® 


-1 


2ai(ai —l)[(rc —l)^^+2m^] t 
^‘^x{x — l)-\-rn^x-\-M'-^ (1~^) -* 


(9.29) 


where C(/x^) = 0 has been considered. On inserting Eqs. (9.27) and (9.29) into Eq. (9.22) and noting that in the 
approximation of order Z 2 — \ should be taken in Eq. (9.22), it can be found that 


= 1 - 


(47i-) 2 3 Jo 


I fo 


= 1 - - - (1 - ^}B 

(ai+1) . 2ai(ai—1) [(ai—1) (/a/m—l)/j.'^+2(/j. —mfi] 

s[td‘^x{x — l)-\-m^x-\-M^{l—x)]^ fi‘^x{x — l)-\-m‘^x-\-M'^{l—x) 


}■ 


(9.30) 


When substituting Eq. (9.30) in Eq. (9.24) and applying the familiar integration formulas, through a lengthy 
calculation, we obtain 


7m(A) — ^^{^i(A) + ^ 2 (A) In ^ + a^^(a)?3(A) 

+ 2^4 K4(A) + x^C5(A)] J(A)} 

where 

? 2 (A) = 3^[fA3 + (6p^-7/32)A2 

- p^n 

6 (A) = f A5 - ( 2 p 2 + 3/32)A4 + i(3p4 + 3p2/32 _ 2/34)y3 

+(P^ - /36(P^ - 4/36A^ - - /36(3p'‘ - V/3^ + P'^)X + {p^ - pP f, 

UX) = f A^ - ( 6 p 2 + 7 / 32 )A^ + i(9p4 + llp 2^2 _ 8^4)^3 
+ ll(p2 - /32)(^2 _ 2^2)y2 _ I (p2 _ ^2)2(3^2 _ ^2)^ + 7(^2 _ ^2)3^ 

^5(A) = ^A^ - (2p4 + 3/34)A6 + i(3p6 + V /34 _ 3^6)y5 

+ (p 2 — /32)(3p4 _ p2p2 _ _ l-(p 2 _ ^2^2^0^4 _|_ 5 p 2^2 _ 3^4^y3 

+ 5 / 32 (p 2 _ ;32)3y2 + 1(^2 _ ^^)4(3^2 _ ^ 2 )y _ (^2 _ ^2)5^ 


(9.31) 

(9.32) 

(9.33) 

(9.34) 

(9.35) 

(9.36) 


KiX) = A^ - 2(/?2 + p2)y2 ^ ^p2 _ ^2)2 

7m(A) = ^{a(A) + 6(A) In f + i^6(A) 
+ 213-[6(A) + 77^6(A)] J(A)} 


(9.37) 


(9.31) 


where 


28 



(9.32) 


6(A) = ^[fA3 + (6p2-7/32)A2 

-|(p"-/3")(3p"-/3")A + 3(p2-/32)2]^ 


6(A) = f A5 - (2p2 + 3/32)A4 + i(3p4 + 3p2/32 _ 2(i^)X^ 
+{p^ - P'^)i.P^ - 4/36A^ --pi.P^- /3^)(3p'‘ - V/3^ + /3'‘)A + (p^ - 


6 (A) = f A5 - (6p2 + 7/32)A4 + i(9p4 + llp 2^2 _ 8;34)^3 

+ ll(p2 - /32)(p2 _ 2;32)^2 -ly- _ /?2)^ + 7(^2 _ ^ 2 ) 3 ^ 


e5(A) = ^A^ - (2p4 + 3/34)A® + A(3p6 + V/?4 _ 3^6)^5 
+(p^ — /3^)(3p^ — — 7/3^)A^ — -{p^ — /9^)^(6p^ + 5p‘^P'^ — 3/3"‘)A^ 

+5/32(p" - /3")'A2 + i(p2 _ /3^)4(3^2 _ ^2)^ _ (^2 _ ^2)5^ 


if(A) = A^ - 2(/32 + p2)^2 ^ (^2 _ ^2)2 


7m(A) — ^^{6(A) + 6(A) In I + x^k{x)^3W 
+ 2^4 K4(A) + j^^6(A)] J(A)} 

where 

Cl (A) = + 2^^ ■*" “ /5^]> 


6(A) = ^[fA3 + (6p2-7/36A2 

-|(p"-/3")(3p"-/3")A + 3(p2-/32)2]^ 


C3(A) = ^A5 - (2p2 + 3/32 )a 4 + i(3p4 + 3p2/32 - 2/34)A3 

+ ( p 2 - / 32)(^2 _ 4 ^ 2)^2 _ i ( p 2 _ ^ 2 )( 3^4 _ 4 ^ 2^2 + ^ 4 )^ + (^2 _ ^ 2 ) 3 ^ 


6(A) = f A5 - (6p2 + 7/32)A4 + i(9p4 + llp2^2 _ 8;34)^3 

+ ll(p2 - /32)(p2 _ _ I(p2 _ ^2)2(3^2 _ ^2)^ + 7(^2 _ ^ 2 ) 3 ^ 


C5(A) = ^A^ - (2p4 + 3/34)A6 + i(3p6 + V /34 _ 3^6)^5 

+ (p^ — /3^)(3p^ — p^/3^ — 7/3^)A^ — A(p2 _ ^2^2^0^4 _|_ ^p2p2 _ 3^4^^3 

+ 5 / 36 p " - / 3 ")' A 2 + i ( p 2 _ / 3 ^) 4 ( 3^2 _ ^ 2 )^ _ (^2 _ ^ 2 ) 5 ^ 


(9.33) 

(9.34) 

(9.35) 

(9.36) 

(9.37) 

(9.31) 

(9.32) 

(9.33) 

(9.34) 

(9.35) 

(9.36) 


K{X) = A^ - 2(/32 + p2)^2 ^ (^2 _ ^2)2 

T(X\ = ^ In A^-(/3^+P^)- 

yiX) A 2-(/32 + p2) + y7^- 


(9.37) 

(9.38) 
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With the anomalous dimension given above, the equation in Eq. (9.23) can be solved and gives the effective mass 
for a quark as follows 


TOfi(A) = niRe 


(9.39) 


where 



(9.40) 


This integral is not able to be analytically calculated even though the coupling constant in Eq. (9.31) is taken to be 
a constant. In the large momentum limit (A ^ oo), Eq. (9.31) tends to 

lm{\) « ^A. (9.41) 

6TTp 

In this limit, it is seen that in the time-like momentum space, we have 

mii{X) « mRe~ ^ 0. (9.42) 


Graphically, we only show in Fig. (11) the s quark effective mass mR^X) given by the timelike momentum subtraction 
with the coupling constant being taken to be a constant. For other quarks, the behavior of their effective masses is 
similar. In Fig. (II), the solid line represents the effective mass given by the massive QCD, while the dashed line 
represents the effective mass given by the massless QCD. From the figure, we see that the effective mass with a finite 
gluon mass behaves as a constant equal to the mass niR in the region [0,1] of A and then rapidly falls to zero when 
A goes from unit to infinity. The effective mass with zero-gluon mass has a peak around A = 1.38. The appreciable 
difference between the both effective masses occurs in the region [0,10] of A. The effective quark mass given by the 
spacelike momentum subtraction can directly be written out from Eqs. (9.31)-(9.40) by replacing the A in 'jm(X) 
with iX and hence the TOfl(A) becomes a complex function. By numerical calculations, it is found that either the real 
part or the imaginary part of the TOfl(A) behaves as an oscillating function with a damping amplitude as exhibited 
in Fig. (12). In the figure, the solid and the dashed lines represent respectively the real part and imaginary part of 
the effective quark mass which is given by the massive QCD, while the dotted line shows the real part of the effective 
quark mass which is given by the massless QCD. It is noted that in the most of practical applications to the both of 
scattering and bound state problems, only the effective quark mass given by the timelike momentum subtraction is 
concerned . 


X. CONCLUSIONS AND DISCUSSIONS 

In this paper, it has been shown that as the massless QCD, the massive QCD established on the basis of gauge- 
invariance has a set of BRST-transformations under which the effective action and generating functional are invariant. 
From the BRST-invariance, we derived a set of W-T identities satisfied by the generating functionals for full Green 
functions, connected Green functions and proper vertex functions. Furthermore, from the above identities, we derived 
the W-T identities respected by the gluon propagator, the three-line and four-line proper gluon vertices and the 
quark-gluon proper vertex. Based on these identities we discussed the renormalization of the propagators and the 
vertices. In particular, from the renormalized forms of the W-T identities obeyed by propagators and vertices, the 
S-T identity for the renormalization constants is naturally deduced. This identity is helpful for the renormalization 
by means of the renormalization group approach. To show the renormalizability of the massive QCD, the one-loop 
renormalization is performed by the renormalization group method. In this renormalization, the analytical expressions 
of the one-loop effective coupling constant, gluon mass and quark mass have been derived. Since the renormalization 
was carried out by employing the mass-dependent momentum space subtraction scheme and exactly respecting the 
W-T identities, the results obtained are faithful and allow us to discuss the physical behaviors of the effective coupling 
constant and masses in the whole range of momentum (or distance). Particularly, the previous result given in the 
minimal subtraction scheme for massless QCD is naturally recovered in the large momentum limit. 

As shown in sections 8 and 9, in the mass-dependent renormalization, it is necessary to distinguish the results 
given by the timelike momentum subtraction from the corresponding ones obtained by the spacelike momentum 
subtraction. For example, one can see from Figs. (4) and (5) that the effective coupling constants given in the 
timelike and spacelike subtraction schemes have different behaviors in the low and intermediate energy region although 
in the large momentum limit, the difference between the both coupling constants disappears. Obviously, the both 
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results obtained in the timelike and spacelike momentum subtraction schemes are meaningful and suitable for different 
physical processes. For instance, when we study the quark-quark scattering taking place in the t-channel, the transfer 
momentum in the gluon propagator is spacelike. In this case, it is suitable to take the effective coupling constant 
and gluon mass given by the spacelike momentum subtraction. If we investigate the quark-antiquark annihilation 
process which takes place in the s-channel, since the transfer momentum is timelike, the effective coupling constant 
and gluon mass given in the timelike momentum subtraction scheme should be used. It is also seen from sections 8 
and 9 that the gluon mass gives a considerable effect on the behaviors of the effective coupling constant and particle 
masses. In particular, the gluon mass plays an crucial role to determine the singular or analytical behavior of the 
effective coupling constant. Since the gauge-invariance does not exclude the gluon to have a mass, it is interesting 
to examine the gluon mass effect on physical processes. At present, the massless QCD has widely been recognized 
to be the candidate of the strong interaction theory and has been proved to be compatible with the present high 
energy experiments. However, we think, the massive QCD would be more favorable to explain the strong interaction 
phenomenon, particularly, at the low energy region because the massive gluon would make the force range more 
shorter than that caused by the massless gluon. As for the high energy and large momentum transfer phenomena, as 
seen from the massive gluon propagator, the gluon mass gives little influence on the theoretical result in this case so 
that the massive QCD could not conflict with the well-established results gained from the massless QCD in the high 
energy domain. 

At last, we would like to make some remarks on the nilpotency problem of the BRST-external sources. In comparison 
of the BRST-transformations and the W-T identities for the massive QCD with those for the massless QCD, it is seen 
that they formally are almost the same. The only difference is that the BRST-transformation for the ghost particle 
field C°'(x) written in Eq. (2.25) has an extra term proportional to cr^, the ghost particle mass squared. Due to 
this term, there also appear some cr^—dependent terms in the W-T identities as denoted in Eqs. (3.18), (4.2), (4.8), 
(5.11) and (6.18). But, in the physical Landau gauge where cr = 0, all the a^-dependent terms disappear. In this 
case, all the BRST-transformations and W-T identities for the massive QCD are formally identical to those for the 
massless QCD. As one knows, for the massless QCD, the composite field functions as defined in Eq. (3.3) have 
the nilpotency property: = 0 under the BRST-transformations [13-16] which guarantee the BRST-invariance 

of the BRST-source terms introduced in the generating functional. This nilpotency property is still preserved for the 
massive QCD established in the Landau gauge because in the Landau gauge, the BRST-transformations are identical 
to those for the massless QCD. However, for the massive QCD set up in arbitrary gauges, we find (IA$i ^ 0, the 
nilpotency loses due to nonzero of the ghost particle spurious mass a. In this case, as pointed out in section 3, 
to ensure the BRST-invariance of the source terms, we may simply require the sources Ui to satisfy the condition 
denoted in Eq. (3.7). The definition in Eq. (3.7) for the sources is reasonable. Why say so? Firstly, we note that the 
original W-T identity formulated in Eq. (3.2) does not involve the BRST- sources. This identity is suitable to use in 
practical applications. Introduction of the BRST source terms in the generating functional is only for the purpose of 
representing the identity in Eq. (3.2) in a convenient form, namely, to represent the composite field functions in the 
identity in terms of the differentials of the generating functional with respect to the corresponding sources. For this 
purpose, we may start from the generating functional defined in Eq. (3.4) to re-derive the identity in Eq. (3.2). In 
doing this, it is necessary to require the source terms ■u^A'I’i to be BRST-invariant so as to make the derived identity 
coincide with that given in Eq. (3.2). How to ensure the source terms to be BRST-invariant? If the composite field 
functions A<I>i are nilpotent under the BRST-transformation, the BRST-invariance of the source terms is certainly 
guaranteed. Nevertheless, the nilpotency of the functions A$i is not a uniquely necessary condition to ensure the 
BRST- invariance of the source terms, particularly, in the case where the functions A$i are not nilpotent. In the 
latter case, considering that under the BRST- transformations, the functions A<I>i can be, in general, expressed as 
(IA$i = where the are some nonvanishing functions, we may alternatively require the sources Ui to satisfy the 
condition shown in Eq. (3.7) so as to guarantee the source terms to be BRST- invariant. Actually, this is a general 
trick to make the source terms to be BRST-invariant in spite of whether the functions A<I>i are nilpotent or not. As 
mentioned before, the sources themselves have no physical meaning. They are, as a mathematical tool, introduced 
into the generating functional just for performing the differentiations. For this purpose, only a certain algebraic and 
analytical properties of the sources are necessarily required. Particularly, in the differentiations, only the infinitesimal 
property of the sources are concerned. Therefore, the sources defined in Eq. (3.7) are mathematically suitable for the 
purpose of introducing them. The reasonability of the arguments stated above for the source terms is substantiated 
by the correctness of the W-T identities derived in sections 4-7. Even though the identities in Eqs. (4.1) and (4.2) are 
derived from the W-T identity in Eq. (3.8) which is represented in terms of the differentials with respect to the BRST- 
sources, they give rise to correct relations between the propagators and/or vertices. For example, the correctness of 
the relation in Eq. (4.8) can easily be verified by the free propagators written in Eqs. (4.10) and (4 14). These 
propagators are usually derived from the generating functional in Eq. (2.12) by employing the perturbation method 
without concerning the BRST-source terms and the nilpotency of the BRST- transformations. A powerful argument 
of proving the correctness of the way of introducing the BRST-sources is that after completing the differentiations 
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in Eq. (3.8) and setting the BRST-sources to vanish, we immediately obtain the W-T identity in Eq. (3.2) which is 
irrelevant to the BRST-sources. Therefore, all identities or relations derived from the W-T identity in Eq. (3.8) are 
completely the same as those derived from the identity in Eq. (3.2). An important example of showing this point 
will be presented in Appendix where an identity derived from the W-T identities in Eqs. (3.8) can equally be derived 
from the generating functional in Eq. (2.12) which does not involve the BRST-sources. 
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XII. APPENDIX: ALTERNATIVE DERIVATION OF THE W-T IDENTITY 

In this appendix, we give an alternative derivation of a W-T identity without concerning the nilpotency of the 
composite field functions appearing in the BRST-source terms and show that the result is equal to the one obtained 
from . generating functional containing the BRST-external sources. Let us start from the generating functional of 
Green functions given in Eqs. (2.12) and (2.13). For simplicity of statement, we omit the fermion field functions in 
the generating functional and rewrite the functional in the form 

Z[J, e] = X / nA-C.C\exp{iS + * / dS[- (5M“ )2 + A“ 

-b + i f d^xd^yC‘^(x)M^^(x, y)C’’(y)} 

where 

and 

M<^\x,y)=d^,[Vf{x)5\x-y)] 

in which T)^{x) was defined in Eq. (2.9). 

When we make the following translation transformations in Eq. (Al) 

q^{x) C^(x)-f d'^y{M-^Y^{x,y)K^{y) 

CYx) CYx) - J d‘^yKYy){M-^Y^{y,x) 

and complete the integration over the ghost field variables, Eq. (Al) will be expressed as 

Z[J,k,K] = 

where Z[J] is the generating functional without the external sources of ghost fields [13,17] 

Z[J] = ^ j V{A)kF[A]exp{iS + * J d^x[-^{df^A;Y + r^A^]} (A6) 


(A2) 

(A3) 

(A4) 


in which 


Af[A] = detM[A] (A7) 

here the matrix M[A] was defined in Eq. (A3). From Eq. (A5), we may obtain the ghost particle propagator in the 
presence of the external source J 


i/S°'^[x,y,J] = 


6^Z[J,K,K] I 


= i{M Yab[x,y,:;jj]Z[J]. 


(AS) 


The above result allows us to rewrite the W-T identity in Eq. (4.1) in terms of the generating functional Z[J] when 
completing the derivative with respect to k''{y) and setting ^ (x) = Yiu) = Oj 
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where 

D’lfiy) = K"{y) - (AlO) 

Uy 

is the ordinary covariant derivative. On completing the differentiations with respect to the source J, Eq. (A9) reads 

A / V[A]Af [A]exp{iS + ^^d^x[-^ [d^A^f + 

X [/ d^yJ^^^{y)Dl\y){M-^Y-{y, x) - (All) 

= 0 

By making use of Eqs. (A3), (A8) and (AlO), the ghost equation shown in Eq. (4.2) may be written as 

M-[x, A](m-1)^'>[x, 2/, A]Z[J] = rV(x-y)Z[J] (A12) 

When the source J is turned off, we get the equation written in Eq. (2.24) which affirms the fact that the ghost 
particle propagator is just the inverse of the matrix M. 

To confirm the correctness of the identity given in Eq. (All), we derive the identity newly by starting from the 
generating functional written in Eq. (A6) which does not involve the BRST-external sources. Let us make the ordinary 
gauge transformation (5A“ = to the generating functional in Eq. (A6). Considering the gauge-invariance of the 

functional integral, the integration measure and the functional = detM[A], we get [13,17] 

5Z[J] = Y D{A)^f[A]Y d^y[J»^^{y) + M^A'^^iy) 

-^d''Ald^]D'Y{y)eYy)eM^s + i!d^x[-^{di^Aif + (ai3) 

= 0 


According to the well-known procedure, the group parameter 6°'{x) in Eq. (A13) may be determined by the following 
equation [10,14,17] 

M^\x)e\x) = d^{Vf{x)9\x)) = y{x) (A14) 

where A“(x) is an arbitrary function. When setting A“(a;) = 0, Eq. (A14) will be reduced to the constraint condition on 
the gauge group (the ghost equation) which is used to determine the 9°‘{x) as a functional of the vector potential A“(a;). 
However, when the constraint condition is incorporated into the action by the Lagrange undetermined multiplier 
method to give the ghost term in the generating functional, the 9°'{x) should be treated as arbitrary function according 
to the spirit of Lagrange multiplier method. That is why we may use Eq. (A16) to determine the functions 9°'{x) in 
terms of the function A“(a;) . From Eq. (A14), we solve 

r (x) = J d^x(M-i)“^(x - y)X\y) (A15) 


Upon substituting the above expression into Eq. (A13) and then taking derivative of Eq. (A13) with respect to A“(x), 
we obtain 


A / D{A)Af[A] f d^ylJ'>^(y) + M^A»^(y) 

-id^At(y)d^]D^Y(y)(^~')""(9 - exp{*5+ 
= 0 


According to the expression denoted in Eq. (2.4) and the identity jbcdj^cn jg g^^gy |.g ggg 

A'’>^{y)D»;{y)iM-Y^^iy - x) = A^^{y)dl{M-Y^^{y - x) (A17) 

By making use of the relation in Eq. (AlO), the definition in Eq. (A3) and the equation in Eq. (A12), we deduce 


Xd';^Al{y)d^D^Y^y){M-Y‘^Yy - 

= Xd-Al{y)5\x -y)- M^d';Al{y){M-Y’^-{y - x) 


(A18) 


On inserting Eqs. (A17) and (A18) into Eq. (A16), we obtain an identity which is exactly identical to that given in 
Eq. (All) although in the above derivation, we started from the generating functional without containing the ghost 
field functions and the BRST-sources and , therefore, the derivation does not concern the nilpotency of the composite 
field functions appearing in the BRST-source terms. This fact indicates that the W-T identities derived in section 3 
are correct and hence the procedure of introducing the BRST-invariant source terms into the generating functional is 
completely reasonable. 
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XIV. FIGURE CAPTIONS 

Fig. (1): The one-loop gluon self-energy. The solid, wavy and dashed lines represent the free quark, gluon and 
ghost particle propagators respectively. 

Fig. (2): The one-loop ghost particle self-energy. The lines represent the same as in Fig. (1). 

Fig. (3): The one-loop ghost-gluon vertices . The lines mark the same as in Fig. (1). 
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Fig. (4): The one-loop effective coupling constants a/{(A) given by the timelike momentum space subtraction. The 
solid and dashed lines represent the coupling constants given by the massive QCD and the massless QCD, respectively. 
The dotted line denotes the coupling constant given in the minimal subtraction scheme. 

Fig. (5): The one-loop effective coupling constants a/j(A) given by the spacelike momentum space subtraction. 
The solid and dashed lines represent the coupling constants given by the massive QCD and the massless QCD, 
respectively. The dotted line denotes the coupling constant given in the minimal subtraction scheme. The dashed- 
dotted line represents the one-loop effective coupling constants aFi{X) given by the spacelike momentum subtraction 
for which the gluon mass is taken to be a smaller value. 

Fig. (6): The one-loop effective coupling constants a/{(A) given by the timelike momentum space subtraction. The 
solid, dashed and dotted lines represent the coupling constants given by taking the quark ffavor Nf = 2,3 and 4 
respectively. 

Fig. (7): The one-loop effective gluon masses The solid and the dashed lines represent the effective masses 

given in the timelike and spacelike subtractions respectively. 

Fig. (8): The one-loop quark-gluon vertices. The lines represent the same as in Fig. (1). 

Fig. (9): The one-loop quark-ghost particle vertices. The lines represent the same as in Fig. (1). 

Fig. (10): The one-loop quark self-energy. The lines represent the same as in Fig. (1). 

Fig. (11): The one-loop effective quark masses m/{(A) given by the timelike momentum space subtraction. The 
solid and the dashed lines represent the effective masses given by the massive QCD and massless QCD, respectively. 

Fig. (12): The one-loop effective quark masses TOfl(A) given by the spacelike momentum space subtraction. The 
solid and the dashed lines represent respectively the real part and imaginary part of the effective quark mass which 
is given by the massive QCD. The dotted line shows the real part of the effective quark mass which is given by the 
massless QCD. 
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